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Abstract 

Collective phenomena in strongly nonequilibrium systems interacting with elec- 
tromagnetic field are considered. Such systems are described by complicated non- 
linear differential or integro-differential equations. The aim of this review is to 
show that many nonlinear collective phenomena can be successfully treated by a 
recently developed method called the Scale Separation Approach whose name is 
due to the idea of separating different characteristic space-time scales existing in 
nonequilibrium statistical systems. This approach is rather general and can be ap- 
plied to various nonequilibrium physical problems, several of which are discussed 
here. The problems considered not only serve as illustrations of the method but 
are quite important by themselves presenting interesting physical effects, such as 
Collective Liberation of Light, Turbulent Photon Filamentation, Superradiant Spin 
Relaxation, Negative Electric Current, and Magnetic Semiconfinement of Atoms. 
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1 Introduction 



Strongly nonequilibrium cooperative processes that occur in statistical systems interacting 
with electromagnetic field are described by complicated nonlinear differential and integro- 
differential equations. For treating such difficult problems, a general approach has been 
recently developed called the Scale Separation Approach whose basic idea is to present 
the evolution equations in such a form where it could be possible to separate several 
characteristic space-time scales. In many cases, different scales appear rather naturally 
being directly related to the physical properties of the considered system. 

Since the scale separation approach makes the mathematical foundation for the follow- 
ing applications, we start the review with presenting the basic techniques of this approach. 
Then we demonstrate it by applying the method to different physical problems related to 
strongly nonequilibrium processes occurring under the interaction of electromagnetic field 
with matter. The considered examples not only serve as illustrations of the method but 
are of importance as such since they concern interesting and rather nontrivial physical 
effects. For consideration, those effects are chosen that have been first correctly described 
or predicted by the authors. Among these effects, we would like to emphasize, as the most 
interesting, the following: Collective Liberation of Light, Turbulent Photon Filamentation, 
Superradiant Spin Relaxation, Negative Electric Current, and Magnetic Semiconfinement 
of Atoms. 

The content of the report is as follows. In Section 2 the Scale Separation Approach 
is described. This method makes it possible to solve, or to strongly simplify, many 
complicated systems of nonlinear differential equations, including stochastic and partial- 
derivative equations. The mathematical procedure of solving nonlinear differential equa- 
tions in the following applications is based on this approach. The examples we consider 
have mainly to do with the evolution equations describing strongly nonequilibrium sta- 
tistical systems interacting with electromagnetic fields. We concentrate our attention 
on collective phenomena whose existence as such, as well as their properties, are due to 
nonlinear effects. This is why we constantly have to deal with nonlinear equations. 

Resonant interactions of electromagnetic field with radiating systems are usually de- 
scribing by the Maxwell-Bloch equations, in which one often passes to the momentum 
representation by means of Fourier transform. But we prefer to work in the Real-Space 
representation, outlined in Section 3, which seems to be more convenient for employing 
the Scale Separation Approach. Another convenient trick we employ is the elimination of 
electromagnetic field from evolution equations. For this purpose, the operator Maxwell 
equations, supplemented by the Coulomb calibration, can be rewritten in the integral 
form connecting the vector potential with the retarded current formed by the radiating 
system. Substituting this vector potential into evolution equations eliminates from them 
electromagnetic field. In this way, we come to the system of equations not containing ex- 
plicitly electromagnetic field, instead of which there appears an effective dipole interaction 
of radiating atoms. After eliminating electromagnetic field, we have less equations, al- 
though the price for this is that these equations become integro-differential. nevertheless, 
the obtained equations are more convenient for applying to them our method of solution. 
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Another important advantage of the derived equations is the possibility of taking into ac- 
count quantum effects. Such effects are often principal, while the standard semiclassical 
Maxwell-Bloch equations cannot take account of them. To simplify evolution equations, 
not loosing quantum effects, is the idea of the Stochastic Mean-Field Approximation of 
Section 4. Since cooperative electromagnetic phenomena are directly related to arising 
coherence, Section 5 gives the definitions for Dynamical Characteristics of Coherence. 

The equations derived in the previous sections and the method of solution developed 
above are applied to several concrete systems exhibiting interesting physical properties. 
In Section 6, we suggest the theory of Collective Liberation of Light, which can occur in 
materials with polariton band gap. In Section 7, we consider the influence of external 
fields on radiation properties of resonant atoms, checking whether it is feasible to get Am- 
plification by Nonresonant Fields. Section 8 discusses the so-called Mossbauer Magnetic 
A nomaly observed in some magnetic materials. In Section 9 the Problem of Pattern Selec- 
tion is analyzed. This problem arises, for instance, when one needs to describe resonant 
media with spatially nonuniform electromagnetic structures. For treating the problem, 
we have suggested an original approach based on probabilistic analysis of possible spatio- 
temporal patterns. This method is applied, in Section 10, to describing Turbulent Photon 
Filamentation in resonant media. 

Scale Separation Approach, being a general method, can be employed for treating 
strongly nonequilibrium systems of different physical nature. In Section 11, it is used for 
giving a thorough picture of Superradiant Spin Relaxation occurring in nonequilibrium 
nuclear magnets. This method also makes it possible to analyse nonlinear differential 
equations in partial derivatives. Such an analysis helps in finding conditions under which 
unusual nonlinear effects can happen. This is illustrated in Section 12 by describing a 
transient effect of Negative Electric Current in nonuniform semiconductors. Another novel 
effect of Magnetic Semiconfinement of Atoms is described in Section 13. Both these effects 
have been predicted by the authors. In Section 14, we discuss conditions when Nuclear 
Matter Lasing could be possible. 

Throughout the review, we consider several physical systems of rather different na- 
ture, Because of this, it is more appropriate to give all details and to discuss the related 
literature in the corresponding sections, limiting the Introduction by a brief enumeration 
of the considered problems. Section 15 contains Conclusion summarizing main results. 

2 Scale Separation Approach 

Because of the pivotal role of this approach for treating physical problems in the following 
sections, we need to start by presenting its general scheme. It is possible to separate five 
main steps, or parts, of the approach: (i) stochastic quantization of short-range corre- 
lations; (ii) separation of variables onto fast and slow; (iii) averaging method for multi- 
frequency systems; (iv) generalized expansion about guiding centers; and (v) selection of 
scales for space structures. Below, these steps are explicitly explained. 
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2.1. Short— Range Stochastic Quantization. When considering nonequilibrium 

processes in statistical systems, one needs to write evolution equations for some averages 
< Ai > of operators -Aj(t), where t is time and % = 1,2,...,N enumerates particles 
composing the considered system. For simplicity, a discrete index % is employed, although 
everywhere in what follows one could mean an operator A{f^ t) depending on a continuous 
space variable fj. 

There exists the well-known problem in statistical mechanics consisting in the fact 
that writing an evolution equation for < A{ > one does not get a closed set of equations 
but a hierarchical chain of equations connecting correlation functions of higher orders. 
Thus, an equation for < Ai > involves the terms as J2j < AiBj > with double correla- 
tors < AiBj >, and the evolution equations for the latter acquire the terms with triple 
correlators, and so on. The simplest way for making the system of equations closed is by 
resorting to the mean-field type decoupling < A^Bj > — > < Ai >< Bj >. When consid- 
ering radiation processes, this decoupling is called the semiclassical approximation. Then 
the term J2j < AiBj > reduces to < Ai > J2j < Bj >, so that one can say that < Ai > is 
subject to the action of the mean field J2j < B 3 , >. The semiclassical approximation de- 
scribes well coherent processes, when long-range correlations between particles govern the 
evolution of the system, while short-range correlations, due to quantum fluctuations, are 
not important. However, the latter may become of great importance if there are periods 
of time when the long-range correlations are absent. For example, this may happen at 
the beginning of a nonequilibrium process when long-range correlations have had yet no 
time to develop. Then neglecting short-range correlations can lead to principally wrong 
results for the whole dynamics. 

To include the influence of short-range correlations, the semiclassical approximation 
can be modified as follows: 

£ < ^ >=< Ai > ^ < B, > , (1) 

where £ is a random variable describing local short-range correlations. It is natural to 
treat £ as a Gaussian stochastic variable defined by its first, <C £ 3>, and second, <C |£| 2 3>, 
moments. According to the short-range character of local fields, we should set 

« £ » = . (2) 

The second moment, aiming at taking into account incoherent local fluctuations, can be 
defined by means of the following reasoning. Consider the equality 



E < A i B 3 > 



> 



< Ai > 



E < B 3 > 



+ « l£l 2 » 



resulting from definitions (1) and (2). On the other hand, wishing to take into account 
both long-range coherent as well as short-range incoherent terms, one should write 

2/2 \ 



E < > 



< A, > 



E < B j > 



+ EI<^> 
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where the first term in the brackets corresponds to the coherent while the second term, 
to incoherent parts. Comparing the latter two equalities, we come to the conclusion that 



«|£| 2 » = £l<^>l 2 - (3) 

j 

As far as short-range correlations and fluctuations are often due to quantum effects, the 
manner of taking them into account by introducing a stochastic variable £ can be named 
the stochastic quantization. Then the decoupling (1) may be termed the stochastic mean- 
field approximation. A similar kind of approximation has been used for taking account 
of quantum spontaneous emission of atoms in the problem of atomic superradiance [1]. 
Somewhat related ideas have also been used in the stochastic quantization of quantum 
field theory [2]. 

2.2. Classification of Function Variations. Employing the stochastic mean-field 
approximation makes it possible to write down a closed set of stochastic differential equa- 
tions. The next step is to find such a change of variables which results in the possibility of 
separating the functional variables onto fast and slow. Let us consider, first, the variation 
of functions in time. Assume that we come to the set of equations of the form 

du , ds 

in which / = f(e, u, s,£,t), g = g(e, u, s, £, t), and e <C 1 is a small parameter. Equations 
(4) are complimented by initial conditions 

u = uo , s = so (t = 0) . (5) 

Here, for simplicity, we deal with only two functions, u and s, and one small parameter 
e. The whole procedure is straightforwardly applicable to the case of many functions and 
several parameters. 

Let the functions / and g be such that 

lim/^0, limeo = 0. (6) 
Then from Eqs. (4) it follows that 

lim ^ , lim ^ = . (7) 

e-0 dt ' e-0 dt V ' 

This permits us to classify the solution u as fast, compared to the slow solution s. In turn, 
the slow solution s is a quasi-invariant with respect to the fast solution u. Thus, we may 
classify the functions representing the sought solutions onto fastly and slowly varying in 
time. 

In the case of partial differential equations, one has, in addition to time, a space 
variable f. Then the notion of fast and slow functions can be generalized as follows [3]. 
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while 



Let r E V, with V being the measure of the volume V, and let t £ [0, T], where T can be 
infinite. If one has 

[ p Qn | pT 

lin ^ < T7/ -^7 ^ > 7^ > lim<-/ V«dt>^0, (8) 

| /■ ^5 | pT 

lim < — / — df > = , lim < - / Vs dt > = , (9) 

£^0 V Jv dt e^O T Jo v ' 

then the solution u can be called fast on average with respect to both space and time, 
as compared to s that is slow on average. In such a case, s is again a quasi-invariant 
with respect to u. In general, it may, of course, happen that one of the solutions is fast 
in time but slow in space, or vice versa, as compared to another solution. Note that 
in the Hamiltonian mechanics quasi-invariant s with respect to time are called adiabatic 
invariants [4] . A generalization of this notion to the case of both space and time variables 
[3] is given by definition (9). 

2.3. Multifrequency Averaging Technique. Let us continue considering the 
ordinary differential equations (4). The generalization to the case of partial differential 
equations can be done similarly to the way discussed at the end of the previous section. 
After classifying the function u as fast and s as slow, we may resort to the Krylov- 
Bogolubov averaging technique [5] extended to multifrequency systems. 

Since the slow solution s is a quasi-invariant for the fast variable u, one considers the 
equation for the fast function, with the slow one kept fixed, 

8X 

— = f(e,X,z,U), (10) 

here s = z being treated as a constant parameter. The initial conditions for Eq. (10) is 

X = u (t = 0). (11) 

The pair of solutions 

X = X{e,z,£,t) , z = const (12) 

are called the generating solutions. Substituting the solution X into the right-hand side 
of the equation for the slow function s, one defines the average 

g(e, z) = < - f g(e, X(e, z, £, t), z, £, t) dt > , (13) 

T JO 

in which r is the characteristic time of fast oscillations. In many cases, it is sufficient to 
set r — > oo. In this way, we come to the equation 

dz 

-^ = eg(e,z), (14) 

with the initial condition 

z = s (t = 0) . (15) 
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The solution to Eq. (14), 

z = z{e,t), (16) 

is to be substituted into X yielding 

y{e,£,t)=X{e,z{e,t\£,t). (17) 

Generating solutions (12) are the first crude approximations one starts with. More elab- 
orate solutions (16) and (17) are termed guiding centers. 

Notice two points that difference the considered way of obtaining the guiding centers 
(16) and (17) from the standard averaging method [5]. The first point is in retaining in 
Eq. (10) the small parameter e, which makes it possible to correctly take into account 
important physical effects, such as attenuation. The standard manner of defining the 
generating solutions with setting e = would result in essentially more rough approxi- 
mations. The second difference is in the occurrence of the stochastic average in definition 
(13), since here we are dealing with stochastic differential equations. 

2.4. Generalized Asymptotic Expansion. The generating solutions (12) play 
the role of the trial zero-order approximation, while the guiding centers (16) and (17) 
essentially improve the trial approximations. Higher-order corrections may be obtained 
by presenting the general solutions as asymptotic expansions about the guiding centers. 
Then, /c-order approximations are written as 

k k 
u k = y(e, f , t) + ]T y n (e, £, t) e n , s k = z(e, f) + £ z n (e, £, t) e n . (18) 

n=l n=l 

Such series are named generalized asymptotic expansions [6], since the expansion coeffi- 
cients depend themselves on parameter e. The right-hand sides of Eqs. (4) are also to be 
expanded about the guiding centers yielding 

k 

f(e, u k , s k , £, t) ~ f(e, y, z, + £ /n(e, f , t) e n , 

n=l 
k 

g(e, u k , s k , £, t) ~ g(e, y, z, f , t) + ^ g n {e, t) e n . (19) 

n=l 

Then, expansions (18) and (19) are to be substituted in Eqs. (4) with equating the like 
terms with respect to the explicit powers of e. Thus, in the first order, this gives 



dyi £ f _ ^ ^ _ / _ N v ( _ j_ ^ dz\ 

where 



d j- = fi(e,t,t) -g(e,z) Xi(e,^,i) , -± = g(e,y, z,£,t) -g(e, z) , (20) 



d 

Xi(e,^,t) = — X(e,z,^t) , z = z(e,t). 
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For the approximations of order n > 2, we get 



dy r 
dt 



■n 



fn(e,C,t) , 



dz r 

dt 



■n 



(21) 



The initial conditions for all n — 1, 2, . . . are 



z n = 



(t 



0). 



(22) 



The functions f n and g n depend on y 1 , y 2 , . . . , y n , and on zi, z 2 , ■ ■ ■ , z n , but it is 
important that the dependence on y n and z n is linear. The latter follows from the fact 
that expanding a function 



and so on. In this way, Eqs. (20) directly define yi and z 1 , and Eqs. (21) are linear 
equations, thus, being easily integrated. 

Usually, one does not need the higher-order approximations since the main physics, in 
the majority of cases, is already well described by the guiding centers (16) and (17). The 
latter are good approximations to the exact solutions [7] in the time interval < t < T s /e, 
where T s is a characteristic time of the slow-solution variation. In those cases when 
the higher-order approximations are important, each /c-order approximant can also be 
improved by invoking some sort of summation [8] of asymptotic series (18), for instance, 
the self-similar summation [9-12]. 

2.5. Selection of Space Structures. The solutions to differential or integro- 
differential nonlinear equations in partial derivatives are generally nonuniform in space 
exhibiting the formation of different spatial structures. And it often happens that a given 
set of equations possesses several solutions corresponding to different spatial patterns or 
to different scales of such patterns [13]. When there is a family of solutions describing 
several possible patterns, the question arises which of these solutions, and respectively 
patterns, is preferable and in what sense could it be preferable. This problem of pattern 
selection is a general and very important problem constantly arising in considering spatial 
structures. In this subsection we delineate a simple way that in many cases helps to solve 
the problem of pattern selection. A more refined theory will be presented in sections 9 
and 10. 

Assume that the obtained solutions describe spatial structures that can be parametrized 
by a multiparameter (3, so that the fc-order approximations Uk{/3,r,t) and Sk(f3,r,t) in- 
clude the dependence on (3 whose value is, however, yet undefined. To define f3, and 




.71 



in powers of e, one has 



h = f'(y)vi , h = \ [f"(y)yi + /'(y)y 2 ] , 



h = h \f"(y)yi + V'\y)v2 + f\y)y 3 ] , 
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respectively the related pattern, one may proceed in the spirit of the self-similar approx- 
imation theory [14-23], by treating (3 as a control function. According to the theory 
[14-23], control functions are to be defined from fixed-point conditions for an approx- 
imation cascade constructed for an observable quantity. For the latter, one may take 
the average energy defined as follows. The internal energy, which is a statistical average 
of the system Hamiltonian, is a functional E[u, s] of the solutions. Taking the A;-order 
approximations for the latter and averaging the internal energy over the period of fast 
oscillations and over stochastic variables, one gets the average energy 

E k {(3) = < - f E[u k (f3, f, t), s k (/3, f, t)] dt > . (23) 

T JO 

For the sequence of approximations, {E k (f3)}, it is possible to construct an approximation 
cascade whose fixed point can be given by the condition 

^ EM = , (24) 

from which one gets the control function /3 — f3 k defining the corresponding pattern. 
According to optimal control theory, control functions are defined so that to minimize a 
cost functional. The latter, in our case, is naturally represented by the average energy 
(23). Hence, when the fixed-point equation (24) has several solutions, one may select of 
them that one which minimizes the cost functional (23), so that 

E k (p k ) =absmin E k {(3) . (25) 

Equations (24) and (25) have a simple physical interpretation as the minimum conditions 
for the average energy (23). However, one should keep in mind that there is no, in general, 
such a principle of minimal energy for nonequilibrium systems [13]. Therefore the usage 
of the ideas from the self-similar approximation theory [14-23] provides a justification for 
employing conditions (24) and (25) for nonequilibrium processes. 

The scale separation approach presented in this section makes it possible to solve 
rather complicated sets of nonlinear differential equations describing various nonequilib- 
rium phenomena in statistical systems. More details on this approach can be found in 
Refs. [24-28]. 

3 Real Space Representation 

When considering the interaction of atoms with electromagnetic fields, one usually em- 
ploys the so-called mode representation, expanding field operators over mode wave func- 
tions [29,30]. These can be either free-mode functions, that is plane waves, or resonator- 
mode functions depending on the resonator geometry. We prefer to deal with the real- 
space representation because of the following reasons: First, the evolution equations in 
this representation are written in a form more convenient for analysing temporal non- 
stationary behaviour of solutions. Second, it is more suitable for describing nonuniform 
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solutions corresponding to self-organized space structures. And third, this representation 
is more appropriate for using the scale separation approach. Since the real space rep- 
resentation is rarely considered in literature, it is worth recalling in brief the derivation 
of the main equations in this representation [31]. To understand the basis of the main 
evolution equations is very important, for these equations will be constantly used in what 
follows. One more peculiarity of the consideration below, differencing it form the standard 
texts, is the comparison of the formulas for the cases of electrodipole and magnitodipole 
transitions. 

Let us have a system of radiators that can be atoms, molecules, nuclei, etc. Assume 
that the size of a radiator, ao, is small as compared to the mean distance between them, 
a, as well as to the characteristic radiation wavelength A, 



^ « 1 , ^ « 1 , 
a A 



(26) 



while the relation between a and A can be arbitrary. Canonical variables related to the 
electromagnetic field are the electric field E and the vector potential A, whose commuta- 
tion relations are 

~E a (r,t), A p {r',t)] = 4ni c 5 a(3 S(f-f'), 

'E a {f,t), E p (f',t)] = [A a {f,t), A p (f',t)] = 0, 

where c is the light velocity and the index a, (5 — 1, 2, 3, or x, y, z, enumerate the Cartesian 
coordinates. The magnetic field is 

H(f,t) = V x A(f,t) . 

To uniquely define the latter, we invoke the Coulomb gauge condition 

V • A(f, t) = . 

Here and in what follows the system of units is used where K = l. 

The radiator charges are described by the annihilation, ip, and creation, field 
operators with the commutation relations 



rl>(r,t), ip f (r',t) =6(f-f'), 



0, 



0, 



V(f,0, E(f',t)\ = [^(f,t), A(f',t) 

in which the indices minus or plus mean the commutators or anticommutators, respec- 
tively, depending on the Bose or Fermi statistics of the charges. 

Assume that in addition to the quantum radiation fields E and H there are classical 
fields Eq and Hq for which we have 

Eo(r, t) = -Vip (r, t) , H (f, t) = V x A (f, t) , V • A (f, t) = . 

These additional fields can be due to external sources or can be created by the matter 
which the radiators are inserted in. 
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Each radiator is also subject to the action of a scalar potential tpi(r) representing all 
stationary Coulomb interactions. Thus, we may introduce the total scalar and vector 
potentials 



N 



<Pt<a(r, t) = ip (r, t) + J2 <Pi(r) , Aotif, t) = A (f, t) + A(f, t) , 



(27) 



i=i 



where N is the number of radiators. Then the local energy operator is defined as 



H(f,t) 



2m 



+ e <ptot(r,t) , 



(28) 



where mo is mass and e, charge of a particle. Omitting here the relativistic term e 2 A 2 ot /c 2 
and using the Coulomb calibration, we have 



H(f,t) 



V 2 ie 

H Aot(r, t) • V + e <ptot(r, t) . 



2m mo c 



The Hamiltonian of the system of radiators interacting with electromagnetic field and 
with matter is written as the sum 



H = H r + Hf + H r f + H m + H, t 



rf 



(29) 



in which the terms represent, respectively, the Hamiltonians of radiators, field, radiator- 
field interaction, matter, and matter-field interaction. The Hamiltonian of the system of 
radiators is 

V 2 N 1 

2m ^ 

This includes also the direct interaction of radiators with matter by means of the effective 
scalar potentials y?j(r). The field Hamiltonian writes 



H r (t) = j^(f,t) 



H 



E 2 (r,t) + H 2 (r,t) 



df 



The radiator-field interaction is described by 



H rf (t)= yV(f,f) 



ie 



m c 



A to t(r,t) ■ V + e<p (r,t) 



ip(f, t) df . 



(31) 



(32) 



The Hamiltonians of matter and of matter-field interaction are to be specified according 
to particular cases under consideration. 

The size of a radiator, according to inequalities (26), is the smallest characteristic 
length. If f is the center-of-mass of a radiator, we shall use the notation 

Ei{t) = E(fi,t) , Hi{t) = H(f h t) , 

Mt) = A(f h t) , E 0i (t) ee E (f h t) , H 0i (t) = H (f h t) . 
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For f in the vicinity of fi, we may write 

<Po(r,t) ~ -f • (f « fi) , 

A (r, H 0i (t) , A(f, t) ~ - I (r - r,) x . 

The energy levels of each radiator are defined by the Schrodinger equation 

V 2 



2mt 



+ e (fi(r) 



i>n(r~ fi) = E n 1p n (f- fi) , 



where it is assumed that all radiators are identical and (pi(f) = y?(r — fi). The eigenfunc- 
tions ipnif — fi) form a complete orthonormal set enumerated by the indices n and i, so 
that 

/ ^mif-^i) ^n{r-fj) df= 5 mn Sij , 5Z^£(r-fi) tl> n (r' -Ti) = S(f-f) . 

in 

With these functions, we may define the density of transition current 



ze 



2m 



and the transition current 



jmn J jmnix) dv . 

We also introduce the electric transition dipole 

d mn = e J ^* m (f) fip n (f) df 
and the magnetic transition dipole 



fimn = ^ j ^ X 3mn{r) df 



(33) 
(34) 

(35) 
(36) 



Using the equality 



V = m 



V 2 

f, - - — + e ipi(r) 
2m 



one can connect the electric transition current (34) and transition dipole (35) as 

— * — * 

jmn i ^mn d mn , ^mn = E m E n . (37) 

The field operators can be expanded over the basis of the wave functions as 

N 

1>(r,t) = J2 J2 C ni(t) 1>n(f- f) ■ 
n i=l 
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From the commutation relations for the field operators one has 

The fact that each radiator is certainly in one of the states labelled by the index n is 
expressed by the unipolarity condition 

Ei(t)U*) = i- (38) 

n 

The wave functions ip n (r — r*j), in agreement with inequalities (26), are localized in a 
small region of the size of a radiator. Such functions are called the localized orbitals. The 
localization condition can be represented by the equality 

J il)*Jf- fi) f(r) ip n (r- fj) df=0 (i ^ j) , 

in which f(f) is a finite function. 

Using the notations and conditions introduced above, we transform the radiator Hamil- 
tonian (30) to the form 

N 

H r (t) = E E E n cl(t) c ni (t) . (39) 

n i=l 

The radiator-field Hamiltonian (32) becomes 



1 



rnn i=\ 

where 



dmn • Eoi(t) H — j mn ■ Ai{t) + f2 mn ■ Bi(t) 



c 



(40) 



Bi(t) = H 0i (t) + Hi(t) (41) 

is the total magnetic field. 

From definitions (34) to (36), we have 

— * — * — * — * —t 

a mn ~ a nm j 3mn ~ Jnm , A* ran ~ f^nm ■ 

Because the wave functions are usually either symmetric or antisymmetric with respect 
to the spatial inversion, so that 

IM-Vl = \M^)\ > (42) 

— * — # 

then we see that d nn = j nn = but, in general, jl nn ^ 0. 

The next approximation that is usually involved is related to the situation when only 
a couple of radiator levels takes part in the considered process. This happens when the 
transition frequency 

uj = uj 2 i = E 2 - E 1 > (43) 

for these two levels is selected by means of an external alternating field whose frequency is 
close to the transition frequency (43). In this way, considering only two levels is equivalent 
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to the quasiresonance approximation. Then, it is convenient to introduce the transition 
operators 

*f (*) = 4(0 c 2i (t) , a+(t) = 4(0 Cli (t) 
and the population-difference operator 



of(t) = 4(t) c 2i (t) - cUt) c u (t), 



so that 



2cl(t) c u (t) = 1 - <7?(*) , 24 (t) C^f) = 1 + <7?(t) . 
The commutation relations for the introduced operators are 

[a~, a+] = -Sij of , [a~, aj\ = af] = , [a~, a*} = 2 5^ a~ , 

[a~, AA = [a~, E 3 ] = [of, A,] = [of, E 3 ] = , 



-2 of , 



where all operators are taken at coinciding times. 
With the notation 

— * — * 

d 2 i = d, fl 21 = fl , 
we have dyi = d*, flu = and consequently 

ji2 = - iuj d* , J 2 i =iuj d. 



(44) 



(45) 



Since only the difference between level energies is measurable, one can set E\ — 0. 
Then the radiator Hamiltonian (39) reduces to 



z i=l 



(46) 



Everywhere in what follows we assume that electromagnetic fields acting on a radiator do 
not change the classification of its energy levels. In the other case it would be impossible 
to talk about quasiresonance. This implies that the interaction energies of a radiator with 
fields are assumed to be much smaller than u . Because of the latter, the term 



N 



o J2 [(#n + #22) + (#22 - #11) of (01 • Bi(t) , 



i=l 



entering the radiator-field Hamiltonian (40), can be neglected as compared to Eq. (46). 
As a result, we obtain 



N 



Hrf(t) = ~ £ 



1 



ji(t) ■ Mt) + di(t) ■ E 0i (t) + jU(t) ■ Bi(t) 



(47) 



where the notation 

ji(t) = i u da+(t) - d* <j~(t) 



di(t) = dat(t)+d* cr-(t) , 
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& (f A** *r(*) (48) 

is used. The Hamiltonian of the matter-field interaction can be written analogously to 
the first term in Eq. (47) as 

i N 

Hmf(t) = -- £ JmAt) ■ Aj(t) , (49) 
° 3=1 

— * 

where N is the number of particles forming the matter and J m j is a local matter current 
having the structure of the operator J m j = (e/m)pj, with pj being the momentum of a 
j-particle. 

The transition between the quantum states tpi and ip2 can be either accompanied by 
the change of parity or not. Then from definitions (35) and (36) it follows that one has 
one of two possibilities: 

d ^ , [1 = (changed parity) ; 

— * 

d — , /t^O (conserved parity) . (50) 

Thus, we actually have to deal with only one of the dipole transitions, either with electric 
or with magnetic. Here we consider them in parallel in order to compare these two cases. 



4 Stochastic Mean— Field Approximation 

Now it is necessary to write down the evolution equations for the operators entering the 
total Hamiltonian (29) whose terms are given by Eqs (46), (31), (47), and (49). The 
Heinserberg equations yield 

\ ^ E(r, t) = V x H(f, t) - ^ J(r, t) , ^ A(f, t) = - E(f, t) , (51) 
which are, actually, the operator Maxwell equations, where the operator of current is 

N N 

Ar,t)=J2 [ji(t)-cf2i(t) x v]<5(r-r0 + £ J m] (t) 5(r - r 3 ) . (52) 
i=i j=i 

For the transition operators we have 

da~ 



dt 



i u 



°i + (&o d-Ai — i d ■ E 0i — i jl ■ B^j of (53) 



for the lowering operator, where ko = ou /c, and the Hermitian conjugate equation for the 
rising operator af. For the population-difference operator we get 

dc z 

-^ = -2k (daf + d* of) • Ai+2 i {daf - d* of) -E 0i +2 i (jl af - /? of) ■B i . (54) 
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From Eqs. (51), using the Coulomb calibration, we find the wave equation 

(*-±£)&,t) = -±A?,t). (55, 

The solution of the latter has the form 

1 r -( If- r'W df' 

A(f, t) = A vac {f, t) + - J J\r',t- ' - ' 1 , (56) 

in which A vac is the vacuum vector potential being a solution of the uniform wave equation. 
With the operator of current (52), the vector potential (56) can be written as the sum 

A = A vac + A rac i + A mat (57) 

— * 

of the vacuum potential A vac , the radiator potential 



A^r„ t) = £ _ I ( ( - ^) + ? | x ( r , _ _ x) ft ( ( _ !k) , (58) 

and of the matter potential 

A mat (f u t) = J2 — Jmj (t~ —) , (59) 

where f^ = fj — r*,-, = \fij\, and the summation J2j does not include the term with 
j = i. 

Our aim is to derive the evolution equations for the variables 

Ui(t) = < a-(t) > , Si (t) = < a*{t) > , (60) 

in which the angle brackets mean the statistical averaging over the radiator degrees of 
freedom. For the double correlators, we shall employ the mean-field type decoupling 

<afaf> = <afxaf> (i ? j) . (61) 

The quantum effects due to self-action [29] can be taken into account by including into 
the evolution equations the attenuation terms defined by 

7 = ^o (^o + ^o) , (62) 

where do = \d\ and /iq = \ More generally, one includes the phenomenological longitu- 
dinal and transverse attenuation parameters 71 and 72. 

To take into account the retardation, we may remember that the action of electro- 
magnetic fields is characterized by the energies that are much smaller than cu . That is, 
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in the zero order one has a { ~ exp(— iu t), as follows from Eq. (53). This suggests to 
treat the retardation by means of the formula 

< crj (t - > = Uj(t) exp(i k r i3 ) , (63) 

which can be called the quasirelativistic approximation since in the relativistic limit c — > 
oo, Eq. (63) becomes an identity. 

Comparing the terms of the vector potential (58), induced by either electrodipole or 
magnetodipole transitions, we notice their essential difference. Really, averaging over 
angles gives 

E fra) r i3 = , (64) 
j 

unless there is a special arrangement of radiators in space. Hence, the vector potential 
induced by magnetodipole transitions, in usual conditions, is negligibly small. Then for 
the averaged potential (58), we have 

< A rad (fi, t) > = iklJ2 (dip* i:j u* - d* tfij Uj) , (65) 

j 

where 

_ exp(i k nj) 
<Pij = 7— • ( 66 ) 

Kq I ij 

The influence of vacuum fluctuations and of matter is characterized by the term 

— * r — * — * -| 

= k d- [A vac (fi, t) + A mat (fi, t)J , (67) 

which we consider as a stochastic variable, whose properties are to be defined by additional 
conditions. 

In this way, we come to the evolution equations for the transverse variable, 
^ = -(i u + 72) Ui - i (d ■ E 0i + \x ■ H 0i ^j + 

+ i kl Si d- E (d^ij u j - d* ifiij Mj) + Si& , (68) 

j 

and for the longitudinal variable, 

^ = 2 i u* (d- E 0i + jl- i? i) - 2iui(d* • E 0i + p* ■ H ^j - 
-2ikl (du* + d* m) ■ E (dtfj u* - d* ^ uj) - 7i (si -C)-2 (u* & + u z C) , (69) 

3 

where ( G [—1, 1] is a pumping parameter. An equation for u* can be obtained by the 
complex conjugation of Eq. (68). Another useful equation is 



d\ui 



dt 



-2 72 \ui\ 2 + Si (u* & + Ui Ci ) - i s t u* (d ■ E 0i + jl ■ H 0i ) + 
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+ i si Ui (d* ■ E 0i + /I* ■ floi) +ik% Si (u* d + u { d* ) • ^ V*j u ) - d* ifij • (70) 

j 

Equations (68) and (70) are basic for describing nonequilibrium collective phenomena in 
radiating systems. The set of assumptions employed for deriving these equations can 
be briefly named the stochastic mean-field approximation since the mean-field type de- 
coupling (61) was used for the radiator correlators, but quantum effects are taken into 
account through the stochastic variable (67). 



5 Dynamical Characteristics of Coherence 

One of the most important results of the cooperative behaviour of radiators is the ap- 
pearance of coherent radiation. The level of coherence of electromagnetic fields can be 
described by the corresponding correlation functions [32]. Here we introduce another char- 
acteristic of coherence, which is convenient for considering the radiation from ensembles 
of radiators [33]. 

The energy density of the radiated electromagnetic field is 



I (* 2 + " 2 ) • < 71 > 



W - 

8tt 

where E = E(f,t) and H = H(f,t). Differentiating Eq. (71) with respect to time, using 
the Maxwell equations (51), and defining the intensity of scattering 



dW„ 



dt 

and the Poynting vector 

— * 

S = 

8tt 



Uj.E + E-J) (72) 



S=^- (ExH-Hx E) , (73) 



we obtain the continuity equation 

(W + W s ) + div,S = . (74) 

The intensity of radiation into the unit solid angle is 

I(n,t) = <:n-S(r,t) :> r 2 , (75) 

where n = f/r, r = |f|, and the colons imply the normal ordering of operators. To 
accomplish the latter, one separates the Hermitian operators into their conjugate parts, 
which, for instance, for the vector potential (58) reads as 

A rad (f, t) = A + {f, t) + A- (f, t) , (76) 

where 



A + (r,t) = J2 



i kn d 



\r — rj\ 



1 + 


% k 


f 


-fj\ 




r — fj 


3 



jl x (r — fj) 



it - - \ f-fj\ 
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Respectively, the electromagnetic positive and negative fields related to Eq. (76) are 

-* 1 0A rad -*, -*_ -* -* t* -* . -*_ 

Erad = 77— = E + E , H rad = V x A rad = H + H . 

c ot 

In the time and space derivatives, we may employ, for differentiating af, the relations 



In the wave zone, where r ^> \ fi\ and \ f— fj\ ~ r — n • fj, (r ^> \ fj\), we have 

A+{r,t)~i^ (d + fIxn)J2 a t( t - r ~l' ffj ) > ( 7? ) 

from where 

E + = -ik A + , H + = nxE + . (78) 

Then in the part of the Poynting vector (73), describing the radiation from the ensemble 
of radiators, one has 

— * c —* —* —* —* 

Srad ~j E ra d X H rad , Hrad U X E rad . 

For the corresponding part of the radiation intensity (75), we get 

or 2 -* / -* \ 2 
Irad{n, t) = — <■ E 2 rad - (n ■ E rad ) :> . (79) 

Averaging the latter over stochastic variables and over fast oscillations yields 

I(n,t)= / < Irad(n,t) > dt , (80) 

2,71 Jo 



the slow variables in the process of integration being kept fixed. For the radiation intensity 
(79), this results in 



N 



I(n,t) = coo 7 E faW < °t(t)aj(t) > , (81) 

ij 

where 

3 2 

fij(n) = 7T~ \nx e\ exp (ikon- 7%) (82) 

and e = d/d or /7/ /x depending on the type of radiation. 

In the radiation intensity (81), we may separate the terms with the coinciding and 
with different indices, so that Y,ij — J2i=j + X^j- This makes it possible to separate the 
radiation intensity into the incoherent and coherent parts, 

7(n, t) = I inc (n, t) + I coh (n, t) , (83) 
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so that the incoherent radiation intensity is 

1 N 

Iinc(n, t) = - uj 7 £ /«(«) t 1 + ( 84 ) 

2 i=i 

while the coherent radiation intensity is 

N 

I C oh(n,t) = uj 7 fa(n) u *(*) u j(*) ■ ( 85 ) 

Here the equality 2afa~ = 1 + erf was used. The total radiation intensity 
integral 

I(t) = J I(n,t) dn(n) = I inc (t) + I coh (t) 
over solid angles. Here the incoherent part is 

1 N 

kneit) = " WO 7 E t 1 + S *W] > 

^ 1=1 

and the coherent part is 

AT 

W*) =w tE fa «?(*)%(*) , (88) 

where 

fij = J f i:j (n) dn(n) , f u = 1 . (89) 
Finally, the level of coherence can be defined [33] by means of the coherence coefficients 

C coh (n, t) ee , CUt) - Y% ■ 00) 

The radiation is mainly incoherent when C co h <C 1 and it is almost purely coherent if 
C co h 3> 1. 

6 Collective Liberation of Light 

A system of initially inverted atoms can, due to photon exchange, become strongly corre- 
lated, as a result emitting a coherent pulse. This effect of self-organization, accompanied 
by a coherent burst, is called the Dicke superradiance [34]. This phenomenon is well 
studied for atoms in vacuum [1,29,30], including different particular cases, such as super- 
radiance in two-component systems [35-37], superradiance from ensembles of three-level 
molecules [1,38], two-photon superradiance [39,40], and so on (see citations in Refs. [41]). 
When radiating atoms or molecules are placed in a solid, they interact with phonons 
[42,43], which can lead to such interesting phenomena as the laser cooling of solids [44,45]. 



is given by the 
(86) 

(87) 
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When an atom is placed in a periodic dielectric structure, in which, due to periodicity, 
a photonic band gap develops, then spontaneous emission with a frequency inside the 
band gap can be rigorously forbidden [46,47]. This kind of matter, where photon band 
gap appears because of the structure periodicity in real space, has been called photonic 
band-gap materials. The photon band gap also appears in natural dense media due 
to photon interactions with optical collective excitations, such as phonons, magnons, or 
excitons [48,49]. One calls this type of the gap the polariton band gap since photons 
coupled with collective excitations of a medium are termed polaritons. 

If a single resonance atom is placed in a medium with a photon band gap, and the 
atomic transition frequency lies inside this gap, then the spontaneous emission is sup- 
pressed, which is named the localization of light [46,47]. This effect is caused by the for- 
mation of a photon-atom bound state [50-52]. When a collection of identical resonance 
atoms is doped into a medium with a photon band gap, so that the atomic transition 
frequency is inside this gap, then the atoms, in principle, can radiate because of the for- 
mation of a photonic impurity band within the photon band gap [50,53-55]. A model 
case of a concentrated sample, whose linear size L is much smaller than the radiation 
wavelength A, has been considered for studying superradiance near a photonic band gap 
[56,57], when the transition frequency almost coincides with the frequency of the upper 
band edge. Here, following Ref. [58], we study the realistic case of a sample with \ <^ L. 

Assume that the localization of light occurs for a single atom with an electric dipole 
transition, so that its population difference is always s = s(0). Considering an ensemble 
of resonance atoms, we resort to Eqs. (68), (69), and (70). For simplicity, we write Ui = u 
and Si = s. Introduce the effective coupling parameters 

_ 37 ^ sin(A: r i3 ) , _ 37 cos(/c r^) 

9 ~ 4 72 t fc o r y ' 9 ~ 4 72 k ry ' ^ 

where 7 = 4/cq In the absence of resonator imposing a selected mode, 

g^g'^^pX 3 , (92) 

where p is the density of resonance atoms. It is convenient to introduce the effective 
frequency and effective attenuation defined, respectively, as 

Q = uj + 72 g' s , r = 72(1-5-5). (93) 

These expressions include the influence of local fields [59] through the coupling parameters 
(91). Since the latter take into account the existence of an ensemble of atoms, we may 
call VL and Y the collective frequency and collective width, respectively. 
With these notations, Eq. (68) reduces to 

ci ii 

_ = _ (iQ + r) u + s i + 72 e d 2 (g + % g') s u* , (94) 
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— * 

where £ = and e d = d/ofo. Equation (69) becomes 



— = -4 72 M 2 - 7i (s - s ) - 2 («* £ + u C)~ 

- 2 72 [(<? + ^ «/) (u* e d f + (g-i g') {u e d *f] , (95) 

where ( = s takes into account that for a single atom the localization of light occurs. 
And for Eq. (70), we have 



d\u 



dt 



-2 T \u\ 2 + s (it* £ + r «) + 72 s [(g + i g') (u* e d f + (g - ig') (u e/f] ■ (96) 
Let us accept the natural inequalities 

72 



^«1 



< i . 



(97) 



And, as always, we keep in mind that the interaction term (67) is small as compared to 
the frequency Q, or that <C £ ^> = 0, which tells that this term is small on average. Then, 
according to Sec. 2, we may classify the solution u as fast while s and|w| 2 as slow. Solving 
Eq. (94), with s being a quasi-invariant, we get 



U + S f e^+ r )*' £(£') dt' 

Jo 



e -(i n+r) t _ 



(98) 



(99) 



«(f) 

Introduce the notation 

R,e /" r 

a = lim — — / < f(t)«(t) > d* , 

r ^oo T T S Jo 

where Re means the real part and which, if £ 3>= 0, takes the form 

a= lim ^ r d f /' e -(in+r)(f-f) «£*(£)£(*') » ^ . 

When £ (£) is a stochastic variable corresponding to a stationary random process, so that 

« eco » = « rc* - o » , 

then the notation (99) becomes 



a = lim 

T^OO 



^ [ dt f e -(^+ r )*' «r(O£(0)» eft', 
r 1 Jo Jo 

Defining a new function 



w = \u\ 2 — a s 2 , 



(100) 

and averaging the right-hand sides of Eqs. (95) and (96) over time and over stochastic 
variables we get 



ds * , d\u\ 2 

- = -4272^-71 (*-C), ^7" 
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-2 r w , 



where 



7i = 71 + 4 7 2 a , C = — s • 

7i 



In what follows, we assume that the quantity (99), describing the intensity of interac- 
tion between atoms and matter, is small, 



Id < 1 



(101) 



To understand the structure of the atom-matter coupling a, we may model the random 
variable £ by the interaction of an atom with an ensemble of oscillators as 



= E 7. (b u e~** + bt e*") , 



where 6^ and 6^ are Bose operators. Then the atom-matter coupling is 

1 + n w 



a 



£7* 



+ 



with n w = <ti blbu 3>. If the coupling a is small, then 7^ « 71, ~ so 5 and c? | ^ | 2 / tit 
dw/dt. Therefore, we obtain the equations 



ds 



-4 72 iu - 71 (s - s ) , 



-2 72 (1 - 5 s) w 



For transient times, when t <C 7f 1 , Eqs. (102) can be solved explicitly, giving 

't- v 



(102) 



-^tanhf^ + i, 



7o 2 



4^ 2 7l 



sech z 



To 



(103) 



where the integration constants 70 = r 1 and to ar e defined by the initial conditions 
u(0) = u and s(0) = sq. For the radiation width 70, we get the equation 



7o 2 



where 



T = 72(1 - gs ) , 
For the delay time, we find 



ro + Vt| (|«o| 2 - a sg) , 

« = a(0) 

7o - r 



(104) 



7o = — 
to 



+ r ° 1 
r = — m 

2 



Introducing the critical coupling 



70 + r 

l-^o) 2 



Ag 2 sl 



+ 



2 ' 
so 



(105) 



(106) 
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we may rewrite the radiation width as 



70 = 2g \s \ 72 V«c - a ■ (107) 

In the case of only one atom, we have to set g — 0. Then Eqs. (102) give 

s = s , w = (\u \ 2 - a si) e" 272 ' (g = 0) , 

which means that the light is localized. But for an ensemble of atoms the radiation 
becomes possible. 

To find out what happens at large times, when t — > oo, we need to analyse the 
stationary solutions of Eqs. (102). There are two pairs of such solutions: 

si = s , wl = (108) 

and 



1 * _ 7i (ffso ~ 1) 

' W 2 ~ AO 

9 4# 2 72 

The stability analysis [58] shows that the fixed point (108) is stable for gso < 1 and 
unstable for gso > 1, when the point (109) becomes stable. When gso < 1, the stationary 
point (108) is a stable node, while that (109) is a saddle point. In the interval 1 < gs < 
1 + 71/872, the fixed point (108) is a saddle point, and that (109) is a stable node. For 
gs > 1 +71/872, the stationary solutions (108) correspond again to a saddle point, while 
the fixed point (109) becomes a stable focus. In the latter case, the pulsing regime of 
radiation is realized, with the asymptotic period between pulses 

Tp= I7i 2 + 8(l-^ )7i7 2 | 1 / 2 ' (U0) 

However, at finite times the radiation pulses are not periodic, so that the characteristic 
time (110) is an approximate period only for t — > 00. 

In this way, when a single atom cannot radiate because of the localization of light, 
an ensemble of atoms can emit coherent radiation, provided that the interaction between 
atoms is sufficiently strong, so that gs > 1. This is why such an effect can be called 
the collective liberation of light. However, this liberation is not complete but only partial 
since s^ > 0. 



7 Amplification by Nonresonant Fields 

An essential enchancement of radiation can occur due to correlations between radiators, 
which results in the emission of a coherent pulse. In order that these correlations could 
be sufficiently strong, it is usually required that the radiation wavelength would be much 
larger than the mean distance between radiators. If the latter is not the case, it is hardly 
probable that the self-organized coherence can develop. How would it be possible to 
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amplify the radiation intensity for a system of radiators whose wavelenght is smaller than 
or comparable with the mean distance between them? This question is of high importance 
for short-wave emission such as x-ray and 7-ray radiation. Coherent transient effects due 
to phase modulation of recoilless 7 radiation have been considered both theoretically and 
experimentally [60-63]. A regenerated signal of gamma echo has been observed [64], which 
is similar to photon echo in optics [65]. In the present section we explore the conditions 
when stationary enchancement of short-wave radiation is feasible, being due to external 
nonresonant fields. Some preliminary results on the problem have been reported [66- 
68], based on simplified models. Here the problem is considered more accurately, using 
the main Eqs. (68) to (70). The latter, in the case of short-wave radiation, when the 
interaction of radiators can be neglected, take the form 

du ■ — * — * 

= -(iu + 72K - isd ■ E 0i , (111) 

^ = 2i(u*d- u t d*) ■ E 0t - 7l ( Si - C) , (112) 

= -2 72 H 2 - is t {u*d- u t d*) ■ E 0i . (113) 

The initial conditions are -Uj(O) = u and Sj(0) = s . 

Assuming, as usual, the existence of small parameters 

J^«l, (114) 
^0 ^0 

we see that Ui has to be classified as a fast solution while s« and | Hi I , clS slow ones. With 
Si being a quasi-invariant, Eq. (Ill) gives 

u At) = e - {tuJ " +r2) 1 \u -i s t d- f E 0i (t) e ^ 0+ ^ T dr . 

Jo 

— * — * 

Let the external field E 0i = E 0i (t) consist of two parts, 

E 0i = E + E x e l ^- Wt) + Ele-^-^ , (115) 

— * 

one being a stationary nonresonant field E , and another part is a pair of plane waves, 
which are in quasiresonance with the transition frequency, 

I A| , x 

J — L <1, A = lo-uj . (116) 

Then the solution of Eq. (Ill) writes 



= _ ,s 1 d 1 E 1 + SiC^Jh ei(S ^_ wt)+ 
ujq -i"f2 A + 272 
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+ I uo + - e* ^ ) e -(^) * . 

w - h2 A + ?7 2 / 

Substituting this into the right-hand side of Eq. (112) and averaging over time as 



(117) 



1 r T 

lim — / f(s, t) dt , 



we come to the equation 



with 



dsj 
dt 



(118) 



\d-E \ 2 + \d-E l 



u 2 + 7o 2 A 2 + 7 2 



2 ' 



7* = 7i + 472 

The solution to Eq. (118) is 

S ,(t) = s e^ + C(l-e-^) . 
Calculating the correlation function 



c* = — c 



(119) 



ut(t)uj(t) = s 2 (t) 



\d-E \ 2 \d-E 1 



+ 



u 2 + 72 2 A 2 + 7 2 



g 2 k-Tij 



where, for simplicity, we set Sj = s, we find the incoherent and coherent radiation inten- 
sities (84) and (85), respectively, as 



I co h(n,t) 



3N 

167T 

3N 2 
' ~8tT 



u 7 \n x q| 2 [1 + s(t)] , 



W07 |n x ed\ s (t) x 



x 



F(k ft) 



\d-E \ 2 

u 2 Q + 7 2 2 



+ F(k n-k) 



\d-Ex\ 2 
AM^f 



(120) 



where n = f/r and the form factor is 



1 



N 



Ak-fi 



N 



N 

E 

i=i 



;m) 



As is seen from expressions (120) and (121), the maxima of coherent radiation occur in 
the directions satisfying the condition 



(k n-k) ■r i = 2-KUi (rii = 0,1,2,...) . 



(122) 



This corresponds either to forward scattering, when all rii = 0, and the periodicity of 
matter is not required, or to the scattering in the Bragg directions, for which the strict 
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space periodicity of radiators is needed. The enhancement of coherent radiation in the 
directions defined by condition (122) is called the Borrmann effect [69,70], which for the 
case of 7-rays is sometimes termed the Kagan-Afanasiev effect [71,72]. 
The total radiation intensities (87) and (88) are 

hnc{t) = X ~ NCUOI [1 + S(t)] , 

I coh (t) = c 7 s\t) ( + l 4r^) , (123) 



ujI + 7 1 A 2 + 7 | 



where the shape factor is 

3 



ip = —f\nx e d \ 2 F(k n-k) dfi(n) . (124) 
8n J 

The value of the latter strongly depends on the shape of the considered sample. Thus, 
for pencil-like or disk-like shapes [29], one has 



up = 



— -A- <c 1 — <gr 1 



where R and L are the radius and length of a cylindrical sample, and A = 2n/k, k = 

— * 

\k\ = co/c. 

Consider the stationary limit t — > 00, keeping in mind the situation typical of Mosbauer 
experiments, when the alternating field is weak, 



\d-E 1 



|2 

< 1 , (125) 



7i72 

and let us set, for simplicity, 

C = -1 (126) 

which means that there is no additional pumping except through the given field (115). 
Then Eq. (119) reduces to 

;hns , (t) = _ 1 + ^ (\±U + \d-E A 



t^oo ^ \ C^o A 2 + 7I 

For the coherence coefficient, defined in Eq. (90), we get 



lim C coh (t) =N^. (127) 
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The role of the nonresonant field E can be characterized by the switching factor [24] 

S(E ,t)= 1{t) (128) 
lim Eo ^o I(t) 
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and its stationary limit 



For our case, we obtain 



S(E ) = lim S(E ,t) 



S(E ) = 1 + 



A 2 



d- En 



d-E 1 



(129) 



(130) 



The switching factors (128) and (129) show how the radiation intensity is amplified when a 

— * — * 

nonresonant field E is switched on, as compared to the situation when E = 0. As is seen 
from expression (130), the amplification can be quite noticeable only if \d- E \ ^> \d- Ei\, 
so that to compensate the smallness of the parameters |A|/wo and 72/^0- 



8 Mossbauer Magnetic Anomaly 

Stationary fields, electric or magnetic, can be due not to external sources but can arise in a 
sample as a result of phase transitions [73,74]. If an ensemble of radiators is incorporated 
into matter exhibiting a phase transition accompanied by the appearance of a constant 
field, the latter may influence some radiation characteristics. An interesting example 
of this kind is given by the gamma radiation of Mossbauer nuclei placed into magnetic 
materials. This example is especially intriguing because of long-standing controversy 
related to its interpretation. 

There exists a number of experiments demonstrating the so-called magnetic anomaly 
of the Mossbauer effect in materials undergoing magnetic phase transition. This anomaly 
consists in an essential increase, up to 50%, of the area under the Mossbauer spectrum 
below the temperature of magnetic transition, as compared to the spectrum area in param- 
agnetic state above the transition temperature. A detailed discussion of these experiments 
can be found in the book [75] and review [76]. The controversy related to this anomaly 
concerns the explanation of the cause of the latter. 

The area of the Mossbauer spectrum, for Mossbauer nuclei in a solid sample, is given 
by the integral 

/+00 
a ahs {uj) du , (131) 
-00 

in which 

f M = exp(-A; 2 r 2 ) (132) 

is the Mossbauer factor, ko = uo/c, r is the mean-square deviation of the nucleus from 
a lattice site, 



(u - u ) 2 + r 2 abs 



<ww) = 7 — ^ (133) 



) 2 + rl 

is the absorption cross-section, T a b s is the absorption half-width, 



2tt(1 + 2h) 

a ° ~ k 2 (l + 2I )(l + a e ) (134) 
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is the cross-section of resonant absorption, Iq and I\ are the nuclear spins of the ground- 
state and excited levels, and a e is the electron conversion coefficient. After integrating 
Eq. (131), we have the spectrum area 

Aabs = KfMVoTabs ■ (135) 

It is important to emphasize that the Mossbauer anomaly, we consider here, has been 
observed only in the so-called absorption geometry, when absorbing Mossbauer nuclei 
are placed inside magnetic matter which is irradiated by an external source. Contrary to 
this, in the experiments with the so-called source geometry, when a radioactive source is 
incorporated into the magnetic matter, but absorbing Mossbauer nuclei are outside this 
matter, no magnetic anomaly has been observed [77-79]. Therefore it is clear that the 
considered Mossbauer anomaly is directly related to the action on Mossbauer nuclei of an 
effective magnetic field appearing below the critical point. But what is the origin of this 
anomaly? 

Historically, the first suggestion was to ascribe the anomaly in the temperature be- 
haviour of the spectrum area (135) to the influence of the appearing magnetic order on 
the Mossbauer factor (132). A number of citations having to do with this suggestion are 
listed in Refs. [75,76]. This assumption implies that the mean-square deviation defin- 
ing the Mossbauer factor (132) is essentially influenced by arising magnetic order. The 
course of reasoning is as follows. Mossbauer nuclei doped into a solid are characterized 
by the same mean-square deviation as the particles forming the solid sample. The latter 
can be described by the Hamiltonian 

Hrn = E f£ + \ E $ (^') - E £ ■ Sj , (136) 

in which $(Rij) is a potential of direct pair interactions while I (Rij) is that of exchange 
interactions, Si is a spin operator, and R^ = \Rij\, with R^ = Ri — Rj. The indices of 
summation in Eq. (136) run as % — 1, 2, . . . , Nq, with Nq being the number of lattice sites. 
Introduce the deviation from a lattice site, 

Ui = Ri- di , (137) 

defined so that 

a { = <Ri> , < Ui > = . (138) 

Taking into account that \ui\ is small as compared to the interparticle distance, one 
expands the interaction potential in powers of uf up to the second order, which results 
in the Hamiltonian 

H m = U + H p + H s + H sp + H' , (139) 
whose terms are explained below: the constant part of the lattice energy 

= \ E ®( a ij) ; a ij = \<*ij\ > = «i - % , (140) 
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the phonon term 



h p = E tt- + 1 E E ®?f < «? > ( 141 ) 

i j^j a/3 



in which $^ = <9 2 $(ajj) / 'dafddj , the spin Hamiltonian 



i7 s = - E 'K) % ; S<j = £ • Sj , (142) 



the term responsible for spin-phonon interactions, 



h sp = -j2 E ' " / ^ • (143) 

i^j a/3 



where = d 2 I{a,ij) / dafda 1 ^ , and the term 



* = H (144) 

related to the striction energy, where the striction force acting on the site % is given by 
the components 

F ? = - ^ E [*M - 2/(oy)5«] • 

The correct definition of the lattice sites in Eq. (138) presupposes that they serve as 
equilibrium positions for particles. This implies that the striction energy is to be zero on 
average, 

< H' > = . (145) 
Then one invokes a kind of the semiclassical approximation 

< ufSij > = <<>< Sij > = , < ufujSij > = < ufuj >< Sij > , 

decoupling the phonon and spin degrees of freedom, which suggests to present the operator 
term in the spin-phonon interaction (143) as 

< u] S tj = < < > + uf u? < > - <ufu?>< S^ > . (146) 
Thus, the matter Hamiltonian (139) can be reduced to 

H m = U + t p + t s , (147) 

with the renormalized terms 

Z7 = f/ + E E Iff <v?v? j ><S ij >, 

i^j a/3 
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i i^j a/3 i^j 



in which the striction energy, because of condition (145), is omitted and the renormalized 
interactions are 

Djf = < - 2lff < S i3 > , J y ee I(a l3 ) +E iff < uf U P > . 

a/3 

The renormalized dynamical matrix Dff defines the effective phonon spectrum uj ks through 
the eigenvalue problem 

m ■ 

where e ks is a polarization vector, the index s labelling polarizations. The spectrum 
and polarization vectors are assumed to be even functions of the wave vector, so that 
Uks = oj_ ks and e ks = e~ks- Polarization vectors form a complete ortho normal basis with 
the properties 

(iks ■ Cks' = $ss> , E e ks e ks — <W ■ 

s 

Expanding the deviation and momentum as 



one transforms the renormalized phonon Hamiltonian to the standard form 

% = H (&L hs + . 

After this, it is straightforward to calculate the correlators 

< U?U P > = i V COth ^ , 

1 3 2N ^ m cu ks 2T 

in which T is temperature. Thus, one gets the mean-square deviation from the equation 

1 V < u a u a > = 1 V 1 Uh 



^ = i E < « > = f -FE 7 coth ^ . (148) 
— bmJ\l ^ io ks 21 



In this way, the influence of magnetic order on the mean-square deviation comes from its 
influence on the phonon spectrum. 

However, the magnitude of the spin-phonon interaction, renormalizing the dynamical 
matrix, is rather small, as compared to the magnitude of direct interactions [76], so that 
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\I°f I ~ 10 3 . Hence, magnetic order cannot influence much phonon frequencies, as 
well as the sound velocity 

c s = hm ^ = - lim V V D?f {k ' e£ s ef s . (149) 

J a/3 

This conclusion is in agreement with all known experiments where phonon characteristics 
have been examined by means of neutron scattering, sound- velocity measurements, elastic 
and thermal investigations. The onset of magnetic order can change the Mossbauer factor 
not more than by 1%, which cannot explain the observed Mossbauer anomaly of the 
spectrum area (135). 

Another explanation was advanced by Babikova et al. [78], supposing that magnetic 
order can influence the electron conversion coefficient a e in the cross-section (134). A 
noticeable decrease of the conversion coefficient could lead to the increase of the cross- 
section (134), and, consequently, to the increase of the spectrum area (135). The decrease 
of the conversion coefficient could be due to the suppression of the conversion channel in 
favour of the 7-radiation channel whose weight could be increased by the enhancement 
of the 7-radiation caused by the arising magnetic order [80]. 

To estimate the influence of an effective magnetic field, appearing in magnets, on the 
radiation intensity of Mossbauer nuclei, we have to consider the switching factor (130) 
that in our case, takes the form 

S{H Q ) = l + 1 l 
lo 

For the characteristic Mossbauer nucleus 57 Fe, we have o;o = 1.44 x 10 4 eV and 72 = 
71 = 0.67 x 10 _8 eV, which can be transformed to the frequency units as c^o ~ 10 19 s _1 
and 72 ~ 71 ~ 10 7 s _1 . The corresponding wavelength is A ~ 10~ 8 cm. Let us take for 
the effective magnetic field H ~ 10 5 G and for the alternating source field Hi ~ 10~ 5 G. 
The transition magnetic dipole fi ~ 0.1/z n , where \i n is the nuclear magneton, hence 
/i ~ 10~ 13 eV/G. This gives HqHq ~ 10 7 s _1 and fi Hi ~ 10 _3 s _1 . From here we obtain 
^iHl/ulHl ~ 10~ 4 , which tells us that the switching factor S(H ) changes too little. 
Therefore, although the arising magnetic order does enhance the radiation of Mossbauer 
nuclei, this enhancement is not sufficient for causing such a drastic increase of the spectrum 
area. 

The last quantity that could be blamed to be responsible for the Mossbauer magnetic 
anomaly is the absorption width T a b s . The latter can be presented as the sum 

r a6s = 7 2 + 7 2 * (150) 

of the homogeneous line width 72 and of the inhomogeneous line width 73. The inhomo- 
geneous width can be due to the variation of local magnetic fields resulting in the random 
shift of the Mossbauer transition frequency [81]. Returning to Section 3, we see that, really, 
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/7- #0 

it -Hi 



an external magnetic field shifts the transition frequency as cu + (/2 2 2 — A*ii) • H . There- 

— * 

fore, the inhomogeneous width can be of order 7I ~ (/I22 — fin) • Hq or 7^ ~ ^qHq. From 
here, 72 ~ 10 7 s _1 , that is, 7^ ~ 72. In this way, the anomalous increase of the Mossbauer 
spectrum area (135) below the magnetic transition temperature can be explained by the 
increase of the absorption width (150) caused by the increasing inhomogeneous width 

7* ~ HqHq. 

9 Problem of Pattern Selection 

Nonequilibrium cooperative phenomena are often described by nonlinear differential or 
integro-differential equations in partial derivatives. The solutions to such equations are in 
many cases nonuniform in space exhibiting the formation of different spatial structures. 
It happens that a given set of equations possesses several solutions corresponding to 
different spatial patterns [13]. In such a case, the question arises which of these solutions, 
and respectively patterns, to prefer? The problem of pattern selection has no general 
solution [13]. A possible way of selecting spatial structures, by minimizing the average 
energy, was delineated in subsection 2.5. Here we advance another method of pattern 
selection. 

Assume that the considered differential equations in partial derivatives can be reduced 
to a d-dimensional system of ordinary equations; the dimensionality d may equal infinity. 
Suppose also that admissible patterns are parametrized by a multiparameter j3. Let the 
state of the dynamical system be defined by the set 

y(t) = {y t (t)= yi ((3,t)\t = l,2,...,d} (151) 

of solutions to the system of differential equations 

ly(t)=v(y,t). (152) 

For different parameters f3 there are different sets (151) corresponding to different spatial 
structures. All admissible values of (5 form a manifold B = {j3}. Each particular value 
of (3 can be considered as a realization of the random variable from the manifold B. The 
classification of the states (151) can be done by defining a probability measure on B. 

To introduce the probability distribution p(f3, t) of patterns at time t, we resort to the 
ideas of statistical mechanics [82], where a probability p can be connected with entropy 
S by the relation p ~ e~ s . The entropy at time t may be expressed as 

S(t) = ln|A$(£)| (153) 

through the elementary phase volume 

A^(t) = H5 yi (t) . (154) 
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Let us count the entropy from its initial value S(0), thus, considering the entropy variation 

AS(t) = S(t) - S(0) . (155) 
Then the probability distribution p ~ e~ AS , normalized by the condition 

jp(P,t) d(3=l 

takes the form 



e -AS(p,t) 

p(P,t) = — -— , (156) 



Z(t) ' 



where the normalization factor is 



:« = /. 



-AS(J3,t) 



The entropy variation (155) writes 

AS(t) = In 



A<f>(t) 



A$(0) 

where the dependence on /3, for brevity, is omitted. Define the multiplier matrix [83] 

Syi(t) 



(157) 



M(t) = [My (0] , M l3 {t) 



%(0) ' 



for which at the initial time one has 

My(0) = 

The variation of the state (151) gives 



Sy s (0) lj 



Sy(t) = M(t) Sy(0) , 
which yields for the elementary phase volume (154) 

A*(f) = II E My® 8yM- 



(158) 



(159) 



(160) 



* 3 



Hence, the entropy variation (157) is 

AS(t) = In 
With condition (159), this results in 



n e mm 



i 3 



AS(t) = In 



n Mu(t) 



E ^\Mii(t)\ 



;iei) 
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Taking the variational derivative of equation (152), we get the equation 



j t M(t) = J(y,t)M(t) (162) 
for the multiplier matrix (158), where 

J(y, t) = [J^y, t)} , J y (y, t) = , (163) 

is the Jacobian matrix. Substituting the entropy variation (161) into Eq. (156), we get 

n L pw)r; 
p(p>*) = — ^) — ' ( ) 

with 

z(*) = y ni^OMr 1 dp. 

i 

Expression (164) defines the probability distribution of patterns labelled by a multipa- 
rameter p. This expression naturally connects the notion of probability and the notion 
of stability. Really, the multipliers are smaller by modulus for more stable solutions and, 
respectively, patterns, for which the probability is higher. 

Another form of the distribution (164) can be derived as follows. Introduce the matrix 

L(t) = [L tJ (t)} , Lij(t) = In \M l0 {t)\ . (165) 

Then the entropy variation (161) becomes 

AS(t) = Tr L{t) . (166) 

Since the trace of a matrix does not depend on its representations, we may perform inter- 
mediate transformations of Eq. (166) using one particular representation and returning 
at the end to the form independent of representations. To this end, let us consider a 
representation when the multiplier matrix is diagonal. Because of Eq. (162) with the 
initial condition (159), the matrix M is diagonal if and only if the Jacobian matrix is also 
diagonal. Then from the evolution equation (162) it follows that 

M^it) =ex P y^J ii {y{t'),t') dt'} . 

L u (t) = f Ai(t') dt' , Ai(t) = Re J u (t) , 
Jo 

Tr L(t) = f A(t') dt' , A(0 = E Ht)- 



Hence 



from where 
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We assume that the state (151) is formed of real functions, so that the velocity field in 
the evolution equation (152) is also real. Then the eigenvalues of the Jacobian matrix 
(163) are either real or, if complex, come in complex conjugate pairs. Therefore 

£ Re J u (y,t) = ]T J tl (y,t) = Tr J(y,t) . 

i i 

For the entropy variation (166) we obtain 

AS(t) = f A{t') dt! , (167) 



where 

A{t)=TrJ(y,t) (168) 

is called [84] the contraction rate. The latter is given by the form independent of repre- 
sentations of the Jacobian matrix (163). With the entropy variance (167), the probability 
distribution (156) becomes 

= ^ exp {- / o ' A(/M0 dt'} , (169) 
where the contraction rate is defined in Eq. (168) and 



Z(t) = J exp{-^A(/M0 dt'} d/3 . 



The most probable pattern at a time t corresponds to the maximum of the distribution 
(169), 

abs max p(f3,t) -> (3{t) . (170) 
One may also define the average pattern at t as corresponding to 

P(t) = J{3p(P,t)df3. 

The most probable and average patterns, in general, do not coincide, although this may 
happen, especially with increasing time. To illustrate the latter, consider a particular case 
when the contraction rate A((3,t) = A(/3) does not depend on time. Then, as t — > oo, we 
have 

« = / <j( "'* a S ,,, ' A( * ) ' } ' 

where (3q is the point of the minimum of A(/3), so that 

A m = o , A»(/3) = ^ m > o ((3 = (3 ). 
In the distribution 



p(P,t)~ J^^te W {-[A((3)-A((3 )]t} 
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one may expand A(/3) near (5 — flo, which gives 



p(P,t) 



1 



cxp 



( 



(g ~ A>) 2 
la 2 (t) 



I 



a(t) 



A"(A>) * 



1 



From here one finds 



lim p(0,t) = 6(0 - fa) . 

t— >oo 



In this way, if differential equations describing a nonequilibrium process have several 
solutions corresponding to different spatial patterns, the latter can be characterized by 
the probability distribution (169), with the contraction rate (168). In the case when 
the multiplier matrix (158) can be calculated, one may use the expression (164) of the 
probability distribution. If all patterns correspond to stable solutions, it is sufficient to 
analyse only the beginning of the process of pattern formation. Then for the entropy 
variation (167) we may write 



Consequently, the most probable pattern, defined by the maximum of the probability 
distribution (169), that is, by the minimum of the entropy variation (167), is now char- 
acterized by the minimum of the contraction rate A(/3,0) at the initial time. 



Spatial structures can appear in radiating systems if the radiation wavelength is much 
shorter than the system characteristic sizes [13]. For instance, electric field in laser cavities 
can exhibit a state which bears some analogy with a superfluid vortex [85] . The Maxwell- 
Bloch equations for slowly varying field amplitudes have been shown to be analogous 
to hydrodynamic equations for compressible viscous fluid [86]. The Fresnel number for 
optical systems plays the role similar to the Reynolds number for fluids. In the same way 
as when increasing the Reynolds number, the fluid becomes turbulent, there can appear 
optical turbulence when increasing the Fresnel number. 

Spatial structures emerge from an initially homogeneous state with a break of space- 
translational symmetry. For small Fresnel numbers F < 5, such structures correspond to 
the empty-cavity Gauss-Laguerre modes imposed by the cavity geometry. These trans- 
verse structures can be described by expanding fields over the modal Gauss-Laguerre 
functions [87-92], which results in reasonable agreement with experiments for C0 2 and 
Na 2 lasers. For large Fresnel numbers F > 10, the appearing structures are very differ- 
ent from those associated with empty-cavity modes. The modal expansion is no longer 
relevant at large F, and the boundary conditions have little or no importance. The laser 
medium looks like divided in a large amount of parallel independently oscillating uncorre- 
lated filaments [93-100] the number of filaments being proportional to F, contrary to the 
case of small Fresnel numbers when the number of bright spots is proportional to F 2 . This 
filamentation was observed in Dye and C0 2 lasers, as well as in other resonance media, 



AS(P,t) ~ A(/3,0) t 



10 Turbulent Photon Filamentation 
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even without resonators [101-105]. The same type of patterns arises in active nonlinear 
media, such as photorefractive Bi 12 Si02o crystal pumped by a laser [106-109]. In the 
latter media there are also two types of pattern formation: for small Fresnel numbers, 
the symmetry is imposed through the boundary, while for large Fresnel numbers, the 
symmetry is imposed by the bulk parameters. In the case of large F, there occurs a kind 
of self-organization with spontaneous spatial symmetry breaking [110]. It is possible to 
easily notice a qualitative transition in the behaviour of photorefractive media as well as 
in that of lasers: In low-F regime there are a few modes of regular arrangement of bright 
spots corresponding to the peaks of the Gauss-Laguerre functions in cylindrical geome- 
try, the number of modes being proportional to F 2 . And in the high-F regime there are 
many modes spatially uncorrected with each other, which is typical for spatiotemporal 
chaos, the number of the chaotic filaments being proportional to F. Short-range spatial 
correlation is characteristic for turbulence, this is why one calls the similar phenomenon 
in optics the optical turbulence. 

The theory of self-organized photon filamentation in high-Fresnel-number resonant 
media was suggested in Refs. [33,111-116], where the consideration was based on simplified 
models and only the stationary regime was analysed. The choice of filament radii was done 
by means of the variational principle, as is described in subsection 2.5. Here we present 
a more general and elaborate theory based on the evolution equations (68) to (70), which 
includes the description of temporal behaviour, and for defining the characteristics of 
filaments we employ the method of pattern selection developed in Sec. 9. 

First, it is convenient to pass in Eqs. (68) to (70) to continuous representation replac- 
ing the sums by integrals according to the rule 



where p(r) is the spatial density of radiators. Wishing to return to the localized rep- 
resentation, one makes the replacement p(r) = ^2^L 1 S(f ' — r*j). In the case when the 
structure of matter is of no importance, it can be treated as uniform on average setting 
p(f) = p = N/V. Cooperative optical phenomena are often considered in this represen- 
tation of uniform medium [117]. Let us stress that the uniformity of matter in no case 
requires the uniformity of fields or polarization. The solutions to Eqs. (68) to (70) can 
correspond to highly nonuniform structures. 
Introduce the notation 



N 




E 



i=i 




due to an external electric field and of the term 



— * — * 

frad(r, t) = k <d- A rad (f, t) > 



(173) 
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responsible for the action of other radiators. Taking into account Eq. (65), we have 



/r«d(r,f) = -- i 7 p / [(p(r-?')u(r',t)-e d 2 <p*(f-f')u*(r',t)} df' , (174) 
where the continuous representation is used, and 

e ik \r\ 4 

Then Eqs. (68) to (70) acquire the form 

nil n s 

— = -(iwb + 72 )« + a/ , ^ = -2(«7 + /*«) - 7 i(s - , 

±l = - 2l2 \u\ 2 + s(u*f + ru). (175) 
Notice that from the latter two equations one has 

d 



- ( s 2 + 4| M | 2 ) = -2 7 is(s - C) - 872 



u\ 



We consider a sample of the cylindrical shape typical of lasers. The seed laser field 
defining the cylinder axis is given by the sum of two running waves, 

E (f, t) = Ei e l(kz -^ + E\ e- l(kz ^ l) , (176) 

which selects a longitudinal mode. The radius, R, and length, L, of the cylinder are such 
that the following inequalities are valid: 

^«1, (177) 

where a is the mean distance between radiators and A, wavelength. There are also the 
standard small parameters 

— < 1 , — < 1 , — < 1 , (178) 

LOq UJq LU 

with A = u — ujq being detuning. 

The solutions to Eqs. (175) are not necessarily uniform in the whole volume V = nR 2 L 
of the sample, but may have noticeable values only inside narrow regions of filamentary 
form, while being almost zero outside these filaments. Consider one such filament, and 
let us surround it by a cylinder of radius b so that the magnitude of solutions is an order 
smaller at the surface of this enveloping cylinder than at its axis. If the profile of a 
filament is close to the Gaussian exp(— r 2 /2r 2 ), with r/ being the filament radius, then 

b = V2 In 10 77 . (179) 



39 



In what follows we assume this relation between the radius b of an enveloping cylinder 
and the radius 77 of a filament. 

Suppose that there are Nf filaments in the volume of the sample, the axis of each 
filament being centered at a point {x n ,y n }, with n = 1,2, ...,Nf. Let us present the 
solutions to Eqs. (175) as expansions over enveloping cylinders, 

N f N f 

u{r, i) = £ Mr, t) 6 n (x, y) e ikz , s(r, t) = £ s n (r, t) & n (x, y) , (180) 

n=l n=l 

where 



G n (x, y) = e(b- yf(x - x n )* + {y- y n )^ 



is a unit-step function. The filaments are located randomly in the cross-section of the 
sample, but so that their enveloping cylinders do not intersect with each other. The 
interaction between filaments is small, which follows from Eq. (174). This is why they do 
not form a regular lattice but are distributed randomly. 

The function <p(r) in Eq. (174) oscillates at the distance A, and the solutions u n and 
s n essentially change in the radial direction in the interval b. Assuming that 

I « 1 , (181) 

we may say that, in the radial direction, the function tp(r) is fastly varying in space, as 
compared to the slow variation of u n and s n . For the latter, we define the averages 



u 



(t) = ^r I U n (f,t)df, S(t) = ^r I S n (f,t)df (182) 



over the corresponding enveloping cylinder of the volume V n = irb 2 L, where in the left- 
hand side of Eq. (182) we, for short, do not write the index n. 

The seed field (176) is needed mainly for selecting a longitudinal mode with cylindrical 
symmetry, but the amplitude of this field is small, so that 



\d- E x \ 

72 



< 1 . (183) 



The excitation of radiators is accomplished by means of pumping characterized by the 
pumping parameter ( in Eqs. (175). 

Defining the effective coupling parameters 



9- 3lP 



472 K Jv„ 

37P r 

472 V n Jv n 



r sm[k \r-r'\-k(z-z')\ , 

L — uf^\ — drdr ' (184) 
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and the collective frequency and width, respectively, 

Q = u + g'~f 2 s , r = 72(1-55), (186) 

for functions (182) we obtain the equations 

du —* — > 

— = -m + r> - isd ■ E ie - lojt , 

at 

j t = -Ag l2 \u\ 2 - 7i(s - C) " 41m (u* d- E x e"***) , (187) 

Mi = _ 2 I>| 2 + 2s Im (u* d- Ei e~ iu,t ) . 

Because of the inequalities (178) and (183), the solution u in Eqs. (187) is fast, while s 
and \u\ 2 are slow in time. Using the scale separation approach, we find 



sd-E 1 



u(t) = u e- [iU+l)t + " * e -«*_ e -l*»+i;t . (i 88 ) 

uj — \ l + u L J 

Introduce the parameter 

Im f T -» -» 

a = lim — - / u*(t) d ■ E 1 e~ luJt dt , (189) 



characterizing the coupling of radiators with the seed field. This, with Eq. (188), gives 

a = P ■ < 190 > 

The latter, according to inequality (183), is small, 

\a\ < 1 . (191) 

Finally, defining the function 

w=|w| 2 -as 2 , (192) 

we obtain the equations 

ds . . efau 

_ = -4572^-7! (s-C) , — = -2 72 (1 -#s) w . (193) 

The behaviour of solutions to Eqs. (193) essentially depends on the values of the 
coupling parameters (184) and (185). To evaluate the latter, we may notice that their 
integrands diminish and fastly oscillate at the distance of the wavelength A. If condition 
(181) holds, we may neglect boundary effects in the integrals (184) and (185) writing 
approximately 

/ f(r- f ') dfdf' = V n [ /(f) dr. 
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Then parameter (184) reduces to 



9 = 



37T7P 
272 Jo 



r dr 



L / 2 sm(k \/r 2 + z 2 — kz) 



-L/2 



k§\/r 2 + z 2 



dz , 



where r is the radial variable. Because of the quasiresonance condition | A| <C ujq, we have 
k ~ k. With the change of the variable x = k{\/r 2 + z 2 — z), we get 



37T7P 



✓•6 >-fcL 

r dr 

Jo Jkr 2 L 



SIM 



dx . 



2j 2 k JO Jkr 2 /L X 

In this expression, one can replace kL — > 00, thus obtaining 



9 = 



3-7T7P 

272 k Jo 



it / /cr s 
2 " IT" 



where the integral sine appears, 

sint 



7T 



dt = - + S i(a;) , 

or 2 



SI 



i(x) = r 

Joe 



sint 



dt 



Introducing the dimensionless quantity 







A;6 2 2tt6 2 



we come to the coupling parameter 



9 = g(P) = 



37r7pL [P 



- - Si(x) 
2 v ; 



472 k 2 Jo 

This can be integrated explicitly by means of the property 



dx . 



which results in 



j Si(x) dx = x Si(x) + cosrr , 
37T7pL 



472 k 2 







- Si(/3) 



1 — cos/3 j 



For the coupling parameter (185), one similarly finds 

3irjpL ft 3 



9 ' = m 

where the integral cosine occurs, 



472 k 2 



/ Ci(x) dx , 
Jo 



Ci(x) = f 

J oc 



cost 



dt 



(194) 



(195) 



(196) 



(197) 
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Integrating 

/ci(x) dx = xCi(x)-sinx, 

we finally get 

371'jpL 
472 k 2 

To better understand the properties of the coupling parameters, we consider two limiting 
cases. When 1 C 1, then 



JW)= [sm/3-/3Ci(/3)] . (198) 



x 3 „. , s , x 2 



Si(x) ~ x , Ci(x) ~ 7^ + In a; — 

18 4 

where 7^ = 0.577216 being the Euler constant. From here 

37T7pL /7T 1 2 \ , 37T7pL 

f ( x ) - 1 — TT o^-o^ ' 9( x )- -. — p r x\\nx\ 



4 72 A; 2 V2 2 J ' * w 4 72 A; 2 

In the opposite case, when x ^> 1, using 

r~i,/ \ 7r cosx sinx . sin a; cos a; 

Si x ~ — , Ci (x) ~ — , 

2 x x l xx A 

we find 

. . 37T7pL / sinx\ ... 3iT'-fpL /cosx 

9(I) ~ 4^ + — ) • 9 W = *5F I" 

These asymptotic expressions help to analyse the dependence of the coupling parameters 
on the variable (194) changing in the interval 

0<P<2F (^F^^pj . (199) 

The stability analysis of Eqs. (193), similarly to that given in Ref. [58], shows that, 
for g( < 1, the solutions tend to the stationary stable point = (, = 0, while for 
g( > 1, the stable fixed point is 

_ 1 * _ 7i (gC ~ 1) 

5 % 2 7 2 

In this way, for all (3 from the interval (199), except the sole case when g( = 1, there 
exists a stable fixed point, that is, almost all solutions are stable, independently of the 
value of (3. Following the method of pattern selection from Sec. 9, we can equip the 
solutions labelled by (3 with the probabilistic weights (169). The most probable, among 
all stable solutions, is that providing the minimum of the initial contraction rate, which 
for this case is 

A(/?,0) = -71-272(1-^0). (200) 
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The minimum of this rate requires that 



For so > 0, one needs the minimum of g, which gives f3 = 4.9. From Eq. (194), one has 
b = 0.88VXL. And the relation (179) yields 

r f = OAlVXL (s > 0) . (202) 

When So < 0, conditions (201) imply the maximum of g, for which f3 = 1.92, b = 0.55V AL, 
and the filament radius is 

t s = 0.26VAL (s < 0) . (203) 

This is practically the same value as found for the filaments radius in Refs. [33,111-115] 
by using the variational principle of subsection 2.5. When the system of radiators is not 
inverted at the initial time and becomes excited by means of a pulse characterized by the 
pumping parameter (, one has to consider the filament radius (203) as corresponding to 
the most probable pattern. The number of filaments can be defined from the normalization 
condition 

i J s(r,t)df=(, (204) 

assuming that the population difference equals +1 inside each filament of radius 77 and 
— 1 outside of the filaments. Then the number of filaments is 

N f = ^(1 + C)(^) 2 • (205) 

The most probable filament radius (203) and the number of filaments (205) are in good 
agreement with the values observed in experiments [93-99,101-105]. The considered phe- 
nomenon of filamentation can be termed turbulent since the filaments are chaotically dis- 
tributed in space and for sufficiently strong pumping, when g( > 1 +71/872, each filament 
is aperiodically flashing in time. The turbulent photon filamentation is a self-organized 
phenomenon due to the bulk properties of interacting radiators. It practically does not 
depend on boundary conditions and exists in both types of lasers, the resonator-cavity 
lasers, such as CO2 and Dye lasers [93-99], as well as in the resonatorless discharge-tube 
lasers, such as lasers on Ne, Tl, Pb, N 2 , and N2" vapors [101-105]. The turbulent fila- 
mentation is also principally nonlinear phenomenon. Thus, in low-Fresnel-number lasers 
(F < 5) the number of light spots is proportional to F 2 . The same dependence of the 
number of coherent rays on F is typical of the initial linearized stage of superfluorescence 
[118]. However, for high-F lasers (F ^> 10) the number of filaments is proportional to F, 
which is in agreement with formula (205) giving Nf ~ F. 
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11 Superradiant Spin Relaxation 



When the initial state of a spin system is strongly nonequilibrium, different kinds of 
spin relaxation can occur. If there are no transverse external fields acting on spins, they 
relax to an equilibrium state by an exponential law with a longitudinal relaxation time T\. 
When the motion of spins is triggered by a transverse magnetic field, the relaxation is again 
exponential but with a transverse relaxation time T 2 that is usually much shorter than T\. 
A rather special relaxation regime arises, if the spin system is coupled to a resonator. This 
can be done by inserting the sample into a coil connected with a resonance electric circuit. 
Because of the action of resonator feedback field, the motion of spins can become highly 
coherent resulting in their ultrafast relaxation during a characteristic collective relaxation 
time much shorter than T2 [119]. This latter type of collective spin relaxation from a 
strongly nonequilibrium state in the presence of coupling with a resonator is the most 
difficult to realize experimentally and to describe theoretically. Experimental difficulties 
have been overcome in a series of observations of this phenomenon for a system of nuclear 
spins inside different paramagnetic materials [120-127]. The collective relaxation time of 
this ultrafast coherent process is inversely proportional to the number of spins, N, and the 
intensity of magnetodipole radiation is proportional to N 2 , in the same way as cooperative 
radiation time and radiation intensity of iV resonant atoms depend on this number in optic 
superradiance [1,29,30,42,45,59]. This is why the process of collective coherent relaxation 
of spins has been called superradiant spin relaxation or, for short, spin superradiance. In 
the case of spin systems, what is usually measured is not the radiation intensity itself, 
which is rather weak, but the power of current induced in the resonant circuit [128]. The 
enhancement of generated pulses by using resonators is, actually, well known in laser optics 
and is important for realizing superradiance of Rydberg atoms [129] and recombination 
superradiance in electron-hole or electron-positron plasmas [130]. Resonators can be 
employed for modifying radiated pulses in optical superradiance [131]. Note also the 
usage of resonators for amplifying the nuclear spin echo signals in magnets [132,133]. 

The appearance of strong correlations between spins is due to the resonator feedback 
field, but not to the photon exchange as it happens for atomic systems. Hence, various 
quantum effects existing in the interaction of electromagnetic field with atoms [32,134-137] 
seem to be absent in the case of spin systems. Therefore it looked natural to try, for the 
theoretical description of relaxation in a spin system coupled with a resonator, to invoke 
the classical Bloch equations complimented by the Kirchhoff equation for the resonant 
electric circuit [1,119,138-140]. However, these equations can provide a description of 
coherent spin relaxation only when the latter is triggered by a coherent pulse, similarly 
to the semiclassical Bloch-Maxwell equations in optics [1,141,142]. The phenomenon of 
the self-organized coherent spin relaxation cannot be described by the Bloch-Kirchhoff 
equations. Then, what initiates spin motion leading to the appearance of purely self- 
organized spin superradiance? This problem of the origin of pure spin superradiance 
was posed by Bloembergen and Pound [119]. They also noticed that the thermal Nyquist 
noise of resonator cannot be a mechanism triggering the motion of spins, since the thermal 
relaxation time is proportional to the number of spins in the sample and, thus, the thermal 
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damping is to be negligibly small for macroscopic samples. Nevertheless, this notice was 
forgotten by the following researchers who assumed that it is just the thermal noise of 
resonator which triggers the spin motion. 

To resolve this controversy and to discover the genuine mechanisms originating the 
spin motion, it was necessary to turn to microscopic models. The system of nuclear spins 
is characterized [143] by the Hamiltonian 

H=\ E Ha-fr £ B-I t , (206) 

i¥=3 i 

in which spins interact through the dipole potential 



n V r 3 

ij 



where /i n is the nuclear magnetic moment, U is a nuclear spin operator, = \ fij\, = 
fi — fj, Hij = fij/rij. The total magnetic field 

B — H e z + H e x 

contains an external magnetic field H and a resonator feedback field H defined by the 
Kirchhoff equation. 

The temporal behaviour of a finite number of spins, with 27 < N < 343, was analysed 
numerically by computer simulations [144-149]. From various cases studied, we present 
here some that give the general qualitative understanding of the whole picture. In Figs. 
1-4, K coh = Pcoh/Pinc is a coherence coefficient, being the ratio of the coherent part of the 
current power P to its incoherent part, and p z is the negative spin polarization. In Figs. 
5-11, C co h = Icoh/hnc is the coherence coefficient of the average magnetodipole radiation 
defined as in Eq. (90), with respect to the total radiation intensity /. The current 
power and radiation intensity are given in dimensionless units and time is measured in 
units of T 2 . In the figure captions, p z (0) and p x (0) mean the corresponding polarization 
components at the initial time, u>o is the Zeeman frequency, u> is the natural frequency of 
the resonant electric circuit and also a frequency of an alternating magnetic field, if any, 
the amplitude of the latter being denoted by ho- The quantity 

(207) 

is the effective coupling parameter, in which r\ is a filling factor; p n , nuclear density; and 
r2 = T^ 1 is a line width. Computer simulations proved that pure spin superradiance does 
exist with no thermal noise involved. 

However, computer simulations can provide only a qualitative picture, as the number 
of spins considered in such simulations is incomparably smaller than what one has in real 
macroscopic samples. Moreover, these simulations give no analytical formulas, making it 
difficult, if possible, to classify all relaxation regimes occurring when varying the numerous 
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parameters of the system. Simplified models [150] can also provide only a qualitative 
understanding. 

An analytical solution of the evolution equations corresponding to the microscopic 
Hamiltonian (206) and a complete analysis of different relaxation regimes of nonequilib- 
rium nuclear magnets coupled with a resonator has been done [25,26,151-158] by employ- 
ing the scale separation approach. The evolution equations are written for the averages 

« = ^ E < <Sf > , 8 = 1 £ < SI > , (208) 

i i 

where S~ = Sf — iSf. Presenting local fluctuating fields through stochastic variables £o 
and £, one comes [25,26] to the evolution equations 

^ = i (uj - £ + iT 2 ) u - i(~f 3 h + f ) s , 

U S 1 1 

j t = g (-Tsh + u* - - ( l3 h + D u- T^s - C) , (209) 

j \u\ 2 = — 2r 2 |w| 2 - i (73/1 + su* + i (73/1 + C) su , 

in which the resonator feedback field, h, in dimensionless units, satisfies the Kirchhoff 
equation 

dh r* d 

— + 2 l3 h + to 2 / h(t') dt' = -2k - (u* + u)+ 73/ , (210) 
at Jo at 

in which / is an electromotive force, 73 is the resonator ringing width, and k = iir]p n n 2 l /h / y3. 
The random local fields are defined as Gaussian stochastic variables with the stochastic 
averages 

«£o» = «l£| 2 » = l1, (211) 

where is the inhomogeneous dipole broadening. Because of the existence of the small 
parameters 

Ti r 2 r» 73 iai . . 

— < 1 , — < 1 , — < 1 , — < 1 , — < 1 , (212) 

LOq UJq LOq LOo LOq 

where A = uo — loq, the functions u and h can be classified as fast while s and |w| 2 as slow. 

Solving Eqs. (209) and (210), it was shown [25,26,151] that the role of the thermal 
Nyquist noise in starting the relaxation process is negligible. But the main cause triggering 
the motion of spins, leading to coherent self-organization, is the action of nonsecular dipole 
interactions. This gives the answer to the question posed by Bloembergen and Pound 
[119]: what is the origin of self-organized coherent relaxation in spin systems? All possible 
regimes of nonlinear spin dynamics have been analysed. When the nonresonant external 
pumping is absent, that is ( > 0, there are seven qualitatively different transient relaxation 
regimes: free induction, collective induction, free relaxation, collective relaxation, weak 
superradiance, pure superradiance, and triggered superradiance [25,26,151,155]. 
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In the presence of pumping, realized e.g. by means of dynamical nuclear polarization 
directing nuclear spins against an external constant magnetic field, one has ( < 0. Then, 
as was shown using phenomenological equations [139], two stationary solutions can appear. 
In our approach, the behaviour of the system is as follows [158]. When ( < 0, three 
dynamical regimes can be observed, depending on the value of ( with respect to the 
pump thresholds 

Ci = -- , C 2 = -- (l + ^r) ■ (213) 
9 9 V 8IV 

Analysing the equations for the slow variables s and w, where 



r 2 

w = \u\ 2 1 s 2 , (214) 



we find two fixed points 



*; = C, wl = 0; 8^ = --, w; = - Tl{ i + 2 9() . (215) 

9 F 2 g 2 

When (i < ( < 0, the first fixed point is a stable node and the second one is a saddle point. 
For C — Ci) both points merge together, being neutrally stable. After the bifurcation at 
C = Ci, hi the region £2 < ( < Ci, the first fixed point looses its stability becoming a 
saddle point while the second fixed point becomes a stable node. Finally, when C < C2, 
the second fixed point turns into a stable focus, and the first one continues to be a 
saddle point. In this way, there are three qualitatively different lasting relaxation regimes 
induced by the pumping [158]. The first one is the incoherent monotonic relaxation to the 
first stationary solution sj, w\. The second regime is the coherent monotonic relaxation 
to the second stationary solution s* 2) wZ,, although the level of coherence may be rather 
low. And the third case is the coherent pulsing relaxation to the second fixed point. 
This unusual regime of pulsing relaxation was observed experimentally [159]. Here we 
present the results of numerical solution of the evolution equations for the slow variables 
s = z(t) and w(t) defined in Eq. (214). Different cases of the pulsing regime are clearly 
demonstrated in Figs. 12 to 18. In the corresponding figure captions we use the notation 
z Q = z(0), w = w(0), and 7 = 71/72- Everywhere in Figs. 12 to 17, the pump parameter 
is ( = —0.5, and in Fig. 18 this parameter is varied. The coupling parameter (207) is 
always g — 10. 

The problem of superradiant spin relaxation can be generalized to the case of nuclei in- 
corporated into a ferromagnetic matrix, where nuclear and electron spins interact through 
hyperline forces. Some model studies of this case have been undertaken [160-162], and 
a general microscopic theory has also been developed [163]. The latter theory makes it 
possible to discover all feasible causes triggering the process of self-organized coherent 
relaxation. The most important such causes are the dipole hyperfine interactions, dipole 
nuclear interactions, and the transverse magnetocrystalline anisotropy. 
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12 Negative Electric Current 



The study of electric processes in semiconductors is important for describing and mod- 
elling semiconductor devices [164]. One of the most difficult problems is the consider- 
ation of strongly nonequilibrium phenomena in essentially nonuniform semiconductors. 
Nonequilibrium and nonuniform distributions of charge carriers can be formed in several 
ways, for instance, by means of external irradiation [165,166]. Transport properties of 
semiconductors with essentially nonuniform distribution of charge carriers can be rather 
specific. For example, in a sample, biased with an external constant voltage, the result- 
ing electric current may turn against the latter displaying the transient effect of negative 
electric current [3,166-168]. 

Transport properties of semiconductors are usually described by the semiclassical drift- 
diffusion equations [164]. In what follows a plane device, of area A and length L is consid- 
ered, which is biased with a constant voltage V . It is convenient to pass to dimensionless 
quantities, measuring the space variable x in units of L, time in units of the transit time 

L 2 

t = — , n = min{\fa\} , 
faV 

where fa is a mobility of the i-type carriers. And the characteristic quantities 

Q° r\ at? 7T._^0 

Po = , VO = SAI^q , Hjq = — , 

■ _ Qo n _ T/ t _ Po 

3o = — , A) = pV , 4o = — , 

are employed for measuring other physical values which are used below. 
The drift-diffusion equations consist of the continuity equations 

+ ^ + -=6, (216) 

for each type of charge carriers, and of the Poisson equation 

&- 4 'SP. PIT) 

for the electric field E(x,t). Here pi(x,t) is a charge density; fa, Di, and Tj are mobility, 
diffusion coefficient, and relaxation time, respectively; ^ is a generation-recombination 
noise [169]. The sample is biased with an external constant voltage, which in our dimen- 
sionless notation implies that 

/ E(x,t)dx = l. (218) 
Jo 

At the initial time, the distribution of charge carriers 

Pi(x,0) = fi(x) (219) 
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is assumed to be nonuniform. 

The total electric current through the semiconductor sample is 



J(t) = f 1 j(x,t) dx , (220) 
Jo 

where the density of current 

i = E(^-A|) ft + if • (221) 
Because of the voltage integral (218), one has 

J^E(x,t)dx = 0. (222) 

It is also possible to show that 

1 r d 

lim < - / — E(x, t) dt > = . (223) 

r->oo T Jo Ox 

This means that the function E can be considered as slow on average in time and in 
space. Then, treating E as a quasi-invariant, one may find the solutions to Eqs. (216) 
and (217) in order to analyse their general space-time behaviour and to find conditions 
when the effect of negative electric current could arise. Such negative current can appear 
only when the initial charge distribution is essentially nonuniform. For example, if this 
initial charge distribution forms a narrow layer located at the point x — a, then the total 
current (220) becomes negative for a transient interval of time in the vicinity of t = 0, if 
one of the following conditions holds true: 

a < - l — (q > —) , or a>- + (q < — -) , (224) 



where 



Q = Y,Qi, Qi= I pi(x,o) dx . 

A JO 



The effect of the negative electric current can be employed for various purposes, as is 
discussed in Refs. [3,168]. For instance, when the initial charge layer is formed by an ion 
beam irradiating the semiconductor sample, the location a corresponds to the ion mean 
free path. In this case, by measuring the negative current J(0), one can define this mean 
free path 



1 

a — - 



2 4nQ 



AO) 

12i piQi 



(225) 



This formula is valid for both positive and negative values of Q. 

Equations (216) and (217) have also been solved numerically [3,168], which confirmed 
the appearance of the negative electric current. Two cases were analysed, with one layer 
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of charge carriers and with two such layers. Here we present the results of calculations for 
the double-layer case. The initial charge distributions (219) are given by the Gaussians 

in which < a,- < 1 and 




The positive charge carriers, with /ii — 1 and Qi — 1, form the left layer centered at 
ai = a, while the negative charge carriers form the layer centered at a 2 = 1 — a. We keep 
in mind the relation D 2 = 3Di for the diffusion coefficients, typical for holes and electrons, 
and we set D x = 10~ 3 . For short, we use the notation rf 1 = t 2 x = 7 and b x = b 2 = b. 
The generation-recombination noise is neglected, which is admissible at the initial stage 
of the process. As the boundary conditions, we accept the absence of diffusion through 
the semiconductor surface, which implies the Neumann boundary condition 

d 

— Pi(x, t) = (a; = 0, x = 1) . 
ox 

In Figs. 19 to 24, we present the total current (220) as well as the electric current through 
the left surface, J(0, t) = j(0, t) and through the right surface, J(l,t) = j(l,t), defined 
by the current (221) at x — or x — 1, respectively. 

13 Magnetic Semiconfinement of Atoms 

Dynamics of neutral atoms in nonuniform magnetic fields concerns problems of current 
experimental and theoretical interest. By means of such fields, atoms can be confined 
inside magnetic traps, which allows to accomplish various experiments with the systems 
of trapped atoms. Recently, Bose-Einstein condensation has been attained in a dilute gas 
of trapped atoms of 87 Rb [170], 7 Li [171], Na [172], and H| [173]. The details on theory and 
experiment can be found in reviews [174-176]. The Bose-Einstein condensate is believed 
to form, at least partially, a coherent state. If it would be possible to construct a device 
emitting a coherent atomic beam, this would be analogous to a laser radiating a coherent 
photon ray. This is why one may call the device, emitting a coherent atomic beam, an 
atom laser [177-184]. An output coupler, coherently extracting condensed atoms form 
a trap, was demonstrated recently [185-187]. But in these demonstrations, the atoms, 
when escaping from a trap, fly out more or less in all directions, with anisotropy formed 
only by the gravitational force. While the very first condition on a laser is that its output 
is highly directional, with the possibility of varying the beam direction [183]. 

A mechanism for creating well-collimated beams of neutral atoms was advanced in 
Refs. [188-192]. This mechanism suggests an output coupler that extracts trapped atoms 
in the form of a directed beam. 



51 



The motion of neutral atoms in magnetic fields can be described by the semiclassical 
equations for the quantum-mechanical average of the real-space coordinate f = {r a }, 
where a — x,y, z, and for the average S = {S a } of the spin operator [193-195]. The first 
equation writes 

d 2 r a -> dB 

171 ~~dlF = // ° S ' ~dr~ + m9a + ^ a ' ^ 226 ^ 

— * 

where m and /x are mass and magnetic moment of an atom; B is a magnetic field; 
g a is a component of the standard gravitational acceleration; and f a is a collision force 
component. The equation for the average spin is 

— * 

h^- = fi SxB . (227) 

The total magnetic field 

B = Bi + B 2 , 

Bi = B[(x e x + y e y + \z e z ) , B 2 = B 2 (h x e x + h y e y ) , (228) 

where \h\ = 1, consists of the quadrupole field Bi, typical of quadrupole magnetic traps, 
and of a transverse field, e.g., of a rotating field [196,197]. In the quadrupole field, A is 
the anisotropy parameter. 

It is convenient to pass to the dimensionless space variable, measuring the components 
of f in units of the characteristic length 

Bo = §[ ■ (229) 
Introduce the characteristic frequencies by the relations 

dh 



2 ^oB 1 HqB 2 
w, = — — , uj 2 = — — , uj = max 
mRn n t 



dt 



(230) 



Also, we define 



9 a j. fa 



., '6, - (231) 



where 7 is a collision rate and £ a can be treated as a random variable with the stochastic 
averages 

< Ut) > = , < S a (t) {pit') > = 2D a 5 af3 5(t - , 

in which D a is a diffusion rate. Then Eq. (226) can be written as the stochastic differential 
equation 

^ = u\ (S x e x + S y e y + \S Z e z + S) + 7 f , (232) 
and Eq. (227) acquires the form 

— * 

dS 

— = uj 2 As, (233) 
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in which the antisymmetric matrix A = [A a p] has the elements 

A a p — —Ap a , A aa = , 

A 12 = Xz , A 23 = x + h x , A 31 =y + h y . 
Assuming the occurrence of the small parameters 



7 



«1, 



IV i 



UJ 2 



«1, 



^ 2 



«1, 



(234) 



we may classify the variables f and ft, as slow, compared to the fast spin variable S. Then 
Eq. (233) can be solved yielding 



S(t)=J2 Oibi(t) exp{ft(t)}, 



(235) 



i=i 



where 



bi(t) = 



OL\ 2 = i^Ct 



(H = 5(0) • 6<(0) , 
(Ai 2 A 23 - «jA 31 ) e x + (A12A31 + "^23) e y + (A 2 2 + a 2 ) e z 

Q=(Al 2 -\a^) 2 + (Al 2 + \a i \ 2 ) (A^ + A^) , 



a 3 = , 



a 2 ee A? 9 + A 2 3 + A 2 



31 > 



JO 



Substituting Eq. (235) into the right-hand side of Eq. (232) and averaging the latter over 
time and over stochastic variables, we obtain 



dt 2 



F + v\ S , 



(236) 



where 



F ee ivfa 3 < e x + 63 + A63 e z > , 



(x + h° x )S x + (y + h y )S y + XzS° z r _(x + h x ) e x + (y + h y ) e y + Xz e z 

d 3 — — TTTT^ ; TTT^ — ^TTTZ , 03 — 



[(X + /l°) 2 + (y + h° y ) 2 + X 2 Z 2 fl 2 ' 3 [( x + ^)2 + ( y + + A 2 2 2]l/2 > 

angle brackets imply time averaging and h° a = h a (0), S® = S a (0). For the rotating 
transverse field, with 



we find 



h x = cos cut , h y = sin cut , 
B ivl[{l + x)S° x + yS° y + XzS° z ] {xe x + ye y + 2X 2 z e z ) 

r = 



(237) 



2[(1 + 2x + x 2 + y 2 + X 2 z 2 )(l +x 2 + y 2 + X 2 z 2 )] 1 / 2 ' 

The motion of atoms, described by Eq. (236), essentially depends on the initial state, 
which, as is known [198,199], can be prepared in an arbitrary way. Suppose that atoms, 
after being laser cooled in a magneto-optical trap [200], are loaded into a magnetic trap 



53 



where they are further cooled by evaporative cooling down to sufficiently low tempera- 
tures, so that there is a portion of atoms with low velocities, which are located close to 
the trap center. If the initial spin condition for these atoms is such that 5° < and 
S y = S z = 0, then the atoms are confined inside the trap moving in an approximately 
harmonic potential. The gradient of the quadrupole field supplies the levitating force 
to support atoms against gravity. The combination of the magnetic field and gravity 
produces a very nearly harmonic confining potential within the trap volume in all three 
dimensions [201]. 

The semiconfining regime of motion [188-192] can be realized by preparing for the 
spin variable nonadiabatic initial conditions 

S° x = S° y = , S° z = S ^ . (238) 

Such conditions can be arranged in several ways. One possibility could be to confine 
atoms in a trap, where all atoms are polarized having their spins in the z direction, as 
e.g. in the trap of Ref. [201], being a quadrupole trap with a bias field along the z 
axis. Then the longitudinal bias field is quickly switched off, and at the same time, a 
transverse field is switched on, which would correspond to the sudden change of potential 
[202]. Another way could be to prepare spin polarized atoms in one trap and quickly 
load them into another trap with the required field configuration. Atoms can be prepared 
practically 100% polarized [203], with the spin-spin relaxation time reaching 100 s [204]. 
The possibility of realizing two ways of transferring atoms from one trap to another, by 
means of sudden transfer as opposed to adiabatic transfer, is discussed in Ref. [205] . The 
third way of organizing the nonadiabatic initial conditions (238) could be by acting on 
the trapped atoms with a short pulse of strong magnetic field, polarizing atomic spins in 
the desired way. 

With the initial conditions (238), the motion of atoms becomes axially restricted from 
one side, depending on the sign of XS. Atoms fly out of the trap predominantly in one 
direction, forming a well-collimated beam [188-192]. This mechanism can be used for 
atom lasers. Another possibility could be to study the dynamics of binary mixtures of 
Bose systems, where the effect of conical stratification [206] can arise. The mixtures of two 
condensates have been realized for rubidium [207] and sodium [208] , and the dynamics of 
two rubidium condensates was observed in Ref. [209]. 

When solving equation (236) for the realistic case of a finite trap, one should take into 
account the trap shape factor, which can be written in the Gaussian form 

(x 2 -\- y 2 z 2 \ 
"-fP T 2 ) ' 

where R and L are the trap radius and length. The relation between the latter can be 
quite different for different traps, starting from almost spherical traps, where R ~ L, 
to needle-shape traps, with R/L ~ 10~ 3 , as for Ioffe-Prit chard magnetic traps [210]. 
Accepting the initial spin conditions (238), and using the notation 

ip(r) 



[(1 + 2x + x 2 + y 2 + X 2 z 2 )(l + x 2 + y 2 + X^z 2 )} 1 / 2 
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from Eq. (236) we obtain 

^=u> 2 1 Qsfzx + 6 x } > ^ = ul(\*Sfz* + 5 z ) , (239) 

where the equation for y, being similar to that for x, is not written down. Note that 
instead of the Gaussian shape factor for the trap, one could opt for 

ip{r) = 1-B(x-R) G(y - R) 6 (\z\ , 

with O(-) being the unit-step function. 

Equations (239) were analysed both analytically and numerically [188-192]. Their 
solutions display the semiconfined regime of motion. Taking into account random pair 
collisions in Eq. (232) shows that atomic collisions do not disturb the semiconfined motion 
provided that temperature T is sufficiently low, satisfying the condition 

« 1 > (240) 

in which p is the density of atoms and a s , their scattering length. The semiconfined regime 
of motion makes it possible to form well-collimated beams on neutral atoms by means of 
only magnetic fields. 



14 Nuclear Matter Lasing 

The natural question that arises after talking about atom lasers is whether there can be 
produced matter waves corresponding to other Bose particles, which could be employed 
for lasing. One such possibility is related to the creation of large number of pions in 
hadronic, nuclear, and heavy-ion collisions. If the density of pions appearing in the 
course of these collisions is sufficiently high, then correlations between pions can result in 
the formation of coherent state and in the feasibility of realizing a pion laser [211]. Pions 
are not the sole type of Bose particles arising in nuclear matter under extreme conditions 
characteristic of fireballs produced in high-energy collisions [212,213]. There are plenty of 
reviews devoted to the state of nuclear matter at extreme conditions, including the region 
of deconfinement transition. Here we cite only some recent of such reviews [214-217]. 

The very first necessary condition that is required for lasing is to be able to generate 
Bose particles with sufficiently high density. Therefore, in order to answer the question 
what kind of Bose particles appearing in nuclear matter under extreme conditions could 
be used for lasing, one has, first of all, to find out what are these Bose particles and 
under what conditions their density is maximal. In this section, we give a very brief 
account of an analysis based on the multichannel model of nuclear matter [217-221]. The 
main idea in constructing this model goes back to the Weinberg approach for describing 
composite particles [222-224], with effective Hamiltonians that are assumed to be a result 
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of the Fock-Tani transformation [225]. Now we shall not plunge into the details of the 
multichannel model, which can be found in Refs. [217,219], but we shall present some 
figures and will formulate the conclusion of an analysis [221] with regard to the most 
probable candidates for nuclear matter lasing. 

When rising temperature or density, nuclear matter exhibits a transition from hadron 
state to quark-gluon state. This transition is often assumed to be a sharp first-order 
transition. Lattice numerical simulations for the quarkless SU (3) gauge model show that 
deconfinement is really a first-order phase transition [226], which is in agreement with 
the multichannel model. Figures 25 to 27 illustrate the behaviour of some thermody- 
namic characteristics, normalized to the corresponding Stefan-Boltzmann limits, for the 
case of the SU (3) gluon-glueball mixture. Figure 28 shows the related glueball channel 
probability. The sharpness on the deconfinement transition essentially depends on the 
interactions between particles or on their radii, when the composite particles are treated 
as bags [227]. 

In the case of realistic nuclear matter, deconfinement is rather a gradual crossover 
but not a genuine phase transition [217]. Then all thermodynamic characteristics change 
continuously, without jumps. This concerns as well the channel probabilities. Thus, in 
Figs. 29, 30 the channel probabilities of nucleons and dibaryons are shown as functions 
of baryon density normalized to the normal baryon density of nuclear matter hob = 0.167 
fm~ 3 . The possible appearance of dibaryons is of special interest since they, being bosons, 
can form a Bose condensate [217,228-230]. 

Summarizing the results of the analysis [221], three types of Bose particles can appear 
in nuclear matter in large quantities: pions, dibaryons, and gluons. The maximum of the 
pion channel probability, reaching w n = 0.6, occurs in the vicinity of the deconfinement 
transition at T ~ 160 MeV and low baryon densities Ub < Uqb- Dibaryons can appear 
mainly at low temperatures T < 20 MeV and relatively high baryon densities tib ~ 10 n 0B , 
where their channel probability Wq ~ 0.7. Large amount of gluons emerges only at high 
temperatures T > 160 MeV. In addition, one should keep in mind that gluons cannot be 
observed as free particles. 

Talking about possible pion, dibaryon, or gluon lasing from nuclear matter, we have 
touched here just one necessary condition, trying to find out when these Bose particles 
can appear in large quantities. To realize such a lasing in reality will, certainly, require to 
solve a number of other problems. But, anyway, to understand the conditions when this 
lasing could be plausible in principle is the necessary first step. 

15 Conclusion 

We have described a general method for treating strongly nonequilibrium processes in 
statistical systems. This method is called the Scale Separation Approach since its basic 
idea is to try to separate different characteristic scales of time and space variables. The 
idea itself is, of course, not new and we have employed some known techniques. What 
is original in our approach is: (i) The combination of several methods and their adjust- 
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merit to the problems of nonequilibrium statistical mechanics, (ii) The generalization of 
the averaging method to stochastic and partial differential equations, (iii) Probabilistic 
solution of the problem of pattern selection. 

The scale separation approach has been shown to be very useful for describing cooper- 
ative phenomena in the interaction of radiation with matter. To emphasize the generality 
of the approach, it is illustrated here by several different physical examples, whose com- 
mon feature is that the related evolution equations are nonlinear differential or integro- 
differential stochastic equations. Such equations, as is known, are difficult to solve. The 
scale separation approach makes it possible to find accurate approximate solutions. The 
accuracy of these solutions has been confirmed by numerical calculations and by compar- 
ison with experiment, when available. Using this approach, several interesting physical 
problems have been solved and new effects are predicted. Among the most interesting 
applications we would like to emphasize the following. 

Collective Liberation of Light happens when en ensemble of resonant atoms is doped 
into a medium with polariton band gap. If the transition frequency of an atom is inside 
this prohibited gap, then atomic spontaneous emission is strongly suppressed, which is 
termed localization of light. Although spontaneous emission of a single atom is prohibited, 
a collective of such atoms can radiate due to their coherent interactions. As a result of 
this coherent radiation, light becomes partially liberated. We have advanced dynamical 
theory of this light liberation for the realistic situation when the radiation wavelength is 
smaller than the linear sizes of the sample (see Sec. 6). 

Mossbauer Magnetic Anomaly has puzzled researches for many years. This anomaly 
consists in a strong increase of the area under the Mossbauer spectrum, below the tem- 
perature of magnetic phase transition. Several explanations of this anomaly have been 
suggested. We have thoroughly analysed this phenomenon and concluded that previously 
suggested mechanisms cannot explain this anomaly but that its origin is rather in the 
increase of inhomogeneous broadening of Mossbauer nuclei, which is due to the arising 
magnetic field (see Sec. 8). 

Turbulent Photon Filamentation in resonant media is an intriguing example of self- 
organization in a strongly nonequilibrium system, whose dynamical theory was absent. 
We have developed such a theory, based on the probabilistic approach to pattern selection, 
and showed that it gives agreement with experiment (see Sec. 10). 

Superradiant Spin Relaxation occurs in a system of spins coherently interacting with 
each other through resonator feedback field. This ultrafast coherent relaxation is similar 
to superradiance in optical systems, because of which the term spin superradiance was 
coined. Contrary to its optical counterpart, the origin of purely self-organized spin super- 
radiance has not been understood for about 40 years, after Bloembergen and Pound posed 
this problem in 1954. We have developed a theory of nonlinear spin dynamics, based on a 
microscopic Hamiltonian, elucidated the origin of pure spin superradiance, and described 
all main regimes of spin relaxation, without pumping as well as in the presence of the 
latter (see Sec. 11). 

Negative Electric Current is a rather unusual effect, when electric current flows against 
an applied voltage. This is a transient effect that can occur in nonuniform semiconductors. 
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We have predicted this effect and suggested its theory (see Sec. 12). 

Magnetic Semiconfinement of Atoms is another effect we predict. This effect can serve 
as a mechanism for creating well-collimated beams of neutral atoms by means of magnetic 
fields. It can be used to form coherent beams of Bose atoms from atom lasers. We have 
presented a theory of this effect (see Sec. 13). 

The possibility of treating nonequilibrium processes in nonlinear systems of quite dif- 
ferent nature has become possible owing to the Scale Separation Approach, which pro- 
vides accurate approximate solutions to complicated systems of differential and integro- 
differential equations. 
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Figure Captions 



Fig.l. Coherence coefficient K co h, current power P, and spin polarization p z as func- 
tions of time for two different coupling parameters defined in Eq. (207), g\ (solid line) 
and 02 (dashed line), with the relation g\j ' g 2 = 10. 

Fig. 2. The same as in Fig. 1 for two different Zeeman frequencies, Uqi (solid line) 
and u 2 (dashed line), related by the ratio w i/^02 = 5. 

Fig. 3. The same functions as in Fig. 1 for different initial polarizations, p 2 i(0) (solid 
line) and p Z 2(0) (dashed line), with the relation p z i/p Z 2(0) = 2. 

Fig. 4. The same functions as in Fig. 1 for different initial transverse polarizations, 
Pxi(0) (solid line) and ^(O) (dashed line), with the relation p x \/p x i{^) = 0.5. 

Fig. 5. Coherence coefficient C co h, radiation intensivity /, and spin polarization p z 
versus time for p z (0) = 0.48 and different parameters: ujq = 200, g = 25 (solid line); 
ujq = 40, g = 25 (dashed line); and uo = 40, g = 2.5 (solid line with crosses). 

Fig. 6. Coherence coefficient C co h, radiation intensivity J, and spin polarization p z as 
functions of time in the case of switched-off resonator-spin coupling (g = 0). The varied 
parameters are: ujq = 200, p x (0) = 0.48 (solid line); uo = 20, p x (0) = 0.48 (dashed line); 
and c^o = 200, p x (0) = 0.20 (solid line with crosses). 

Fig. 7. The same as in Fig. 6 for p x (0) = 0.48 and for different Zeeman frequencies: 
Uq = 1000 (solid line); u = 200 (dashed line); u = 50 (solid line with crosses); and 
ujq = 200 with switched-off dipole interaction (solid line with triangles). 

Fig. 8. The same as in Fig. 6 for p x (0) = 0.48 but in the presence of an alternating 
magnetic field with the frequency uj = u and different amplitudes: hoi (solid line); /i 2 
(dashed line); where h i/h 2 = 10; and h 03 = (solid line with crosses). 

Fig. 9. The same as in Fig. 8 but for p x (0) = —0.48 and different amplitudes of the 
alternating field: hoi (solid line); /i 2 (dashed line); and /103 (solid line with crosses), where 
the amplitude relations are /ioi/^02 = 0.25 and hoi/h 03 = 0.1. 

Fig. 10. The same as in Fig. 8 for a varying relative detuning from the resonance 
5 = {uj — uq)/uq taking the values: 5 = (solid line); 5 = 0.025 (dashed line); and 
5 = 0.25 (solid line with squares). 

Fig.ll. Radiation intensivity /, coherence coefficient C co h, and spin polarization p z 
versus time, in the absence of alternating external fields and with a weak coupling with a 
resonator, g ~ 1. 
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Fig. 12. Phase portrait demonstrating a stable focus for the parameters Zq = —0.5, 
w = 0.001, g = 10, and 7 = 0.1. 

Fig. 13. Pulsing regime of spin relaxation with the parameters zo = —0.1, wq = 10~ 6 
and 7 = 0.01 for the functions: (a) w(t); (b) z(t). 

Fig. 14. The time dependence of the functions: (a) w(t); (b) z(t), for the parameters 
z = —0.5, wo = 0.001, and 7 = 1. 

Fig. 15. Dynamics of slow solutions: (a) w(t); (b) z(t), for the parameters z = 
-0.5, w = 0.01, and 7 = 0.1. 

Fig. 16. Evolution of slow solutions: (a) w(t); (b) z(t), for the parameters Zo = 
0.5, w = 0.01, and 7 = 0.01. 

Fig. 17. Temporal behaviour of the function w(t) for different sets of parameters: 
(a) z = -0.1, w = 10" 6 , 7 = 0.001; (b) z = -0.1, u = 0.001, 7 = 0.01; (c) 
z = -0.5, wo = 10~ 6 , 7 = 0.1; (d) z = -0.5, w = 0.001,7 = 0.01. 

Fig. 18. Function w(t) for z = 0.5, w = 0.5, 7 = 1, and varying pump parameters: 
C = —0.5 (solid line); ( = —0.3 (dashed line). 

Fig. 19. Electric current through the semiconductor surfaces in the case of a = 
0.1, Q2 = —1, 7 = 1 and different mobilities: ^2 = —10 (solid line); /i 2 = —5 (dashed 
line); /i 2 = —3 (short-dashed line), (a) Left-surface current J(0,t); (b) Right-surface 
current J(l,t). 

Fig. 20. Left-surface current J(0, t) (solid line), right-surface current J(l,t) (dashed 
line), and the total current J(t) (short-dashed line) for a = 0.25, Q 2 = —0.1, /i 2 = — 10 
and different relaxation parameters: (a) 7 = 1; (b) 7 = 10; (c) 7 = 25. 

Fig. 21. Total electric current J(t) for a = 0.25, Q2 = —0.1, /i 2 = —10 and varying 
relaxation parameters: 7 = 1 (solid line); 7 = 10 (dashed line); 7 = 25 (short-dashed 
line). 

Fig. 22. Electric current through semiconductor for the parameters a = 0.1, Q2 = 
— 1, 7 = 1 and different mobilities: fj, 2 = —10 (solid line); /x 2 = —5 (dashed line); /x 2 = —3 
(short-dashed line). 

Fig. 23. Electric current J(t) as a function of time for a = 0.1, /i 2 = —3, 7 = 1 
and different initial charges:: Q2 = (solid line); Q2 = —0.25 (dashed line); Q2 = —0.5 
(short-dashed line); Q2 = —0.75 (dotted line); Q2 = — 1 (dashed-dotted line). 
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Fig. 24. Electric current J(t) for Q2 = — 1, [12 = —3, 7 = 1, and different locations of 
initial charge layers: a = 0.05 (solid line); a = 0.1 (dashed line); a = 0.15 (short-dashed 
line); a = 0.2 (dotted line); a = 0.25 (dashed-dotted line). 

Fig. 25. Relative energy density as a function of temperature in MeV for the SU(3) 
gluon-glueball mixture of different glueball radii: (line 1); 0.5 fm (line 2); 0.7 fm (line 
3); 0.8 fm (line 4); 1 fm (line 5). 

Fig. 26. Relative enthalpy for the gluon-glueball mixture as a function of temperature 
reduced to the deconfinement temperature, in the case of the glueball radius 0.82 fm, 
compared with the lattice numerical calculations. 

Fig. 27. Relative specific heat for the gluon-glueball mixture, for the glueball radius 
0.82 fm, function of temperature in MeV. 

Fig. 28. Glueball channel probability versus temperature in MeV for the glueball radii 
as in Fig. 25. 

Fig. 29. Nucleon channel probability as a function of relative baryon density. 
Fig. 30. Dibaryon channel probability versus relative baryon density. 
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